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Abstract 

A  propulsion  configuration  is  defined  for  a  precision  orbital  tracking  vehicle  (POTV) 
that  employs  electrodynamic  forces  to  control  vehicle  attitude  and  position  with  respect 
to  another  spacecraft.  A  pair  of  electrically  powered  thrusters  assist  the  POTV  in  posi¬ 
tion  control.  The  vehicle  can  maintain  continuous  standoff  at  close  range  from  another 
spacecraft  in  any  direction,  has  docking  capability,  and  has  gross  orbital  transfer  capabil¬ 
ity.  These  capabilities  and  the  specific  design  derived  for  the  POTV  have  been  defined  by 
the  requirements  of  a  postulated  space  facility  (ASSET)  that  salvages  structural  aluminum 
from  the  external  fuel  tank  of  the  Space  Shuttle.  Precision  control  of  the  POTV  is  achieved 
by  an  electrical  conductor  configuration  that  makes  current-produced  thrust  continuously 
available  for  independent  control  of  components  of  vehicle  attitude  and  translation.  A 
robust  tracking  controller  is  designed  which  guarantees  stability  and  compensates  for  geo¬ 
magnetic  field  modeling  error.  This  concept  is  applicable  for  vehicles  with  mass  and  power 
capabilities  differing  from  the  particular  design  described  herein  and  has  application  in  the 
areas  of  attitude  control,  satellite  retrieval,  and  free-flying  platforms. 


AN  ELECTROMAGNETICALLY-CONTROLLED  PRECISION  ORBITAL 


TRACKING  VEHICLE 
(POTV) 


I.  Introduction 

The  use  of  the  earth’s  magnetic  field  as  a  source  of  spacecraft  propulsion  was  first 
proposed  in  1965  by  Drell,  Foley,  and  Ruderman  in  their  ground-breaking  paper  “Drag 
and  Propulsion  of  Large  Satellites  in  the  Ionosphere:  An  Alfven  Propulsion  Engine  in 
Space.”  (2).  Their  idea  has  been  used  to  design  electrodynamic  tethers  for  various  uses 
in  earth  orbit.  Tethers  have  many  non-electrodynamic  uses,  including  payload  orbital 
transfer,  orbital  rendezvous,  docking,  artificial  gravity,  and  gravity-gradient  stabilization. 
These  and  many  other  applications  are  detailed  in  NASA’s  Tethers  in  Space  Handbook 
(15).  Electrically  conducting  tethers  can  be  used  for  power  generation  and  storage,  very 
low  frequency  communication,  orbital  braking,  and  most  importantly  for  our  purposes, 
propulsion.  In  this  thesis  we  design  a  spacecraft  which  uses  the  electrodynamic  forces 
arising  from  the  interaction  of  electrical  currents  with  the  geomagnetic  field  to  precisely 
track  other  space  vehicles  at  close  range,  dock  with  other  spacecraft,  and  perform  gross 
orbital  maneuvers. 

There  mc  two  important  differences  between  the  conductors  the  POTV  uses  to  con¬ 
duct  electricity  uid  the  electrodynamic  tethers  extensively  written  about  in  the  literature. 
First,  there  are  multiple  conductors  attached  to  the  POTV  whereas  the  literature  generally 
discusses  electrodynamic  tethers  singly.  Secondly,  the  POTV  conductors  are  supported  by 
“rigid”  towers  whose  dynamic  behavior  is  much  simpler  tham  that  of  flexible  “strings”.  The 
rigidity  is  provided  by  lightweight  structural  webbing  that  connects  all  the  conductors.  The 
POTV  conductors  are  also  much  shorter  than  the  electrodynamic  tethers  commonly  en¬ 
countered  in  the  literature.  POTV  conductors  are  expected  to  be  less  than  one  kilometer  in 
length.  A  twenty  kilometer  long  electrodynamic  flexible  tether  was  supposed  to  have  been 
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deployed  in  August  of  1992  from  the  space  shuttle.  The  Tethered  Satellite  System  (TSS-1) 
experiment  failed  for  reasons  unrelated  to  the  electrical  characteristics  of  the  tether. 

How  do  the  POTV  conductors  use  the  geomagnetic  field  to  produce  thrust?  We  re¬ 
view  the  pertinent  electromagnetic  theory:  The  force  on  a  current  conductor  in  a  magnetic 
field  is  given  by  the  Lorentz  equation  (15:131): 

t 

Fs=ILxB  (1.1) 

where  F_  is  the  force  exerted  on  the  conductor  by  the  magnetic  field,  I  is  the  current  flowing 
through  the  conductor,  X  is  a  vector  with  magnitude  equal  to  the  conductor  length  which 
points  in  the  direction  of  positive  current  flow,  and  B_  is  the  magnetic  field  vector.  For  our 
application,  H  is  the  geomagnetic  field  vector.  The  Lorentz  equation  is  fundamental  to  the 
thesis  work. 

The  POTV  conductors  must  make  electrical  contact  with  the  plasma  in  the  earth’s 
ionosphere  in  order  to  close  the  current  loop.  Three  basic  methods  for  accomplishing  this 
are  presented  by  Penzo  and  Ammann  (15:120-121): 

•  A  passive  large-area  conductor  at  both  conductor  ends. 

•  A  passive  large-area  conductor  at  one  end  and  an  electron  gun  at  the  other  end. 

•  A  plasma-generating  hollow  cathode  (PMG)  at  both  ends. 

The  first  two  methods  are  dependent  upon  the  density  of  the  surrounding  ionospheric 
plasma.  We  use  PMGs  for  the  POTV  because  this  technology  is  the  most  promising  for 
high  current  densities.  The  power  requirements  inherent  in  electrodynamic  propulsion  are 
discussed  in  appendix  C. 

Since  the  POTV  conductors  rely  upon  the  geomagnetic  field  to  produce  thrust,  the 
magnetic  field  strength  produced  by  the  current  flowing  through  any  one  conductor  at 
another  conductor  should  not  be  a  significant  fraction  of  the  geomagnetic  field  strength. 
The  magnetic  field  (scalar)  due  to  a  current  flowing  through  a  wire  is  (6:243): 


where  /i  is  the  magnetic  permeability  of  free  space,  /  is  the  current  flowing  through  the  wire, 
and  p  is  the  radial  distance  from  the  wire.  This  equation  should  be  checked  when  designing 
the  current  capacity  of  the  conductors  and  the  physical  separation  between  conductors. 

An  electrodynamic  propulsion  systems  is  fundamentally  diflFerent  from  a  conventional 
reaction  systems  in  that  it  produces  thrust  which  is  not  dependent  upon  the  expulsion 
of  mass.  Any  thrust  arising  from  the  PMG’s  ejection  of  electrons/ions  is  assumed  to 
be  negligible.  Unfortunately,  a  major  finding  of  this  thesis  is  that  it  is  not  physically 
possible  to  independently  control  all  components  of  vehicle  position  and  attitude  with 
electrodynamic  thrusters  (conductors)  alone.  Thus  we  add  a  single  pair  of  electromagnetic 
plasma  thrusters  to  provide  needed  position  control  capability.  We  choose  high  specific 
impulse  electrically-powered  thrusters  over  conventional  chemical  thrusters  to  minimize 
fuel  consumption.  Jahn  (8:197)  states  that  plasma  thrusters  should  be  able  to  deliver  high 
thrust  in  conjunction  with  high  exhaust  velocities. 

The  POTV  is  capable  of  gross  orbital  maneuvers  as  well  as  close-range  precision 
tracking.  Penzo  and  Ammann  (15:168-169)  show  how  an  electrodynamic  system  like  the 
POTV  can  be  used  to  change  all  of  the  classical  orbital  elements.  Thus  we  will  not  go  into 
detail  about  large  orbital  change.  We  assume  that  the  POTV  has  enough  gross  orbital 
maneuvering  capability  to  move  close  enough  to  a  target  vehicle  for  the  precision  control 
system  to  handle  close-range  target  tracking  and  docking. 

Docking  using  purely  electrodynamic  forces  has  the  advantage  over  reaction  systems 
of  not  imparting  any  momentum  to  the  tuget  vehicle  from  the  exhaust  products  of  the 
tracking  vehicle.  We  neglect  the  possible  effects  that  the  ionospheric  currents  have  upon 
the  target  vehicle.  There  ace  different  approaches  to  docking  in  the  literature.  Two  in¬ 
plane  methods,  V-BAR  and  R-BAR,  are  described  by  Hall  (5).  We  will  investigate  these 
methods  as  well  as  an  out-of-plsme  approach  (Z-BAR)  to  see  which  approach  is  “best”  for 
the  POTV. 

In  this  chapter  we  have  discussed  how  the  POTV  can  use  the  earth’s  geomagnetic 
field  as  a  source  of  propulsion.  Chapter  two  covers  the  derivation  of  the  relative  trans¬ 
lational  and  rotational  equations  of  motion  for  the  POTV.  In  chapter  three  we  describe 
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the  geomagnetic  field  and  express  it  in  terms  of  an  orbital  frame  of  reference.  Chap¬ 
ter  three  also  covers  magnetic  forces  and  moments.  Chapter  four  covers  the  design  of  a 
POTV  closed-loop  controller.  Chapter  five  goes  into  detail  about  tracking  behavior,  dock¬ 
ing  strategies,  and  conductor-sizing.  The  final  chapter  summarizes  important  results  and 
makes  recommendations  for  further  research.  The  appendices  contain  geomagnetic  field 
plots,  a  description  of  forces  not  included  in  the  ariuations  of  motion,  a  controller  design 
program  listing,  and  plots  corresponding  to  dilFerent  docking  approaches. 
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II.  Derivation  of  Equations  of  Motion 


In  this  chapter  we  derive  the  relative  translational  and  rotational  equations  of  motion 
for  the  POTV.  The  equations  of  motion  derived  herein  include  only  forces  and  moments 
due  to  gravity.  Magnetic  forces  and  moments  are  described  in  chapter  three.  We  describe 
the  coordinate  frames  necessary  for  the  derivation  and  the  rotation  matrices  that  relate  che 
coordinate  systems.  We  also  perform  a  steady-state  analysis  of  the  equations  of  motion. 

2. 1  Coordinate  Frames 

2.1.1  Coordinate  Frame  Description  Figure  2.1  shows  the  coordinate  frames  that 
are  used  in  the  derivation  of  the  relative  equations  of  motion  of  the  tracking  vehicle 
(POTV).  A  description  of  each  coordinate  frame  follows: 

•  Geocentric  Equatorial  Frame  (t):  The  origin  of  the  t  frame  is  at  the  center  of  the 
earth.  The  ti  unit  vector  is  defined  to  point  in  the  vernal  equinox  direction.  The  £3 
unit  vector  is  collinear  with  the  earth’s  rotation  axis  and  points  out  of  the  earth’s 
geographic  north  pole.  The  unit  vector  completes  the  right-handed  orthogonal  set. 
An  excellent  description  of  this  “inertial”  frame  is  found  in  Bate  (1:55). 

•  Orbital  Reference  Frame  (c):  The  origin  of  the  c  frame  is  at  a  point  on  a  reference 
orbit.  This  point  usually,  but  not  necessarily,  coincides  with  the  center  of  mass  of  a 
target  vehicle.  Ci  points  radially  outward  along  a  line  connecting  the  center  of  the 
earth  and  the  c  frame  origin.  62  points  along  the  orbital  velocity  vector  of  the  target 
vehicle.  63  points  “northward”  out  of  the  orbital  plane. 

•  Body-Fixed  Reference  Frame  (6):  The  origin  of  the  b  frame  is  at  the  tracking  vehicle’s 
center  of  mass.  The  coordinate  frame  unit  vectors  coincide  with  the  principal  ajces  of 
inertia  of  the  tracking  vehicle.  This  frame  is  referred  to  as  the  body  frame  henceforth. 

2.1.2  Rotation  Matrices  A  rotation  matrix  is  a  special  type  of  matrix  that  allows 
a  vector  defined  in  terms  of  the  unit  vectors  of  a  coordinate  system  to  be  expressed  in 
terms  of  the  unit  vectors  of  another  coordinate  system.  We  will  derive  the  rotation  matrix 
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TV‘  that  allows  a  vector  expressed  in  inertial  cocM-dinates  to  be  expressed  in  terms  of 
orbital  reference  frame  coordinates.  We  say  that  the  matrix  72*'  transforms  the  vector 
from  the  i  coordinate  system  to  the  c  coordinate  system.  This  matrix  is  of  use  in  section 
3.2  where  we  derive  an  expression  for  the  geomagnetic  field  in  terms  of  orbital  reference 
frame  coordinates. 

The  size,  shape,  and  inertial  orientation  of  the  reference  orbit  and  the  location  of 
the  target  vehicle  within  the  reference  orbit  must  be  known  in  order  to  derive  72*'.  This 
information  is  contained  in  a  set  of  six  orbital  elements  defined  by  Bate  (1:58-60): 

•  The  semi-major  axis  a. 

•  The  eccentricity  e. 

•  The  orbital  inclination  angle  t. 

•  The  right  ascension  of  the  ascending  node  tl. 
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•  The  argument  of  periapsis  u. 

•  The  argument  of  latitude  at  epoch  Uq- 

We  assume  that  the  reference  orbit  is  circular  in  order  to  simplify  the  derivation  of  the 
relative  equations  of  motion.  Thus  the  eccentricity  e  equals  zero  and  the  argument  of 
periapsis  u  is  undefined.  The  orbital  radius  tq  (ro  =  a)  specifies  the  size  of  the  circular 
reference  orbit.  The  right  ascension  of  the  ascending  node  ft  and  the  inclination  i  com¬ 
pletely  determine  the  orientation  of  the  orbit  in  inertial  space.  Note  for  future  reference 
that  t  is  zero  and  that  ft  is  undefined  for  an  equatorial  orbit.  The  argument  of  latitude  at 
epoch  auigle  Uo  locates  the  target  vehicle  within  the  reference  orbit  at  initial  time  (epoch) 
to-  We  define  a  new  angle  u  that  locates  the  target  vehicle  within  the  reference  orbit  at 
any  later  time  t.  The  concept  of  orbital  mean  motion  n  is  used  to  write: 

u(t)  =  «o  +  n(t  -  to)  (2.1) 

where 

where  G  is  the  universal  gravitational  constant  and  is  the  matss  of  the  earth.  The 
angular  velocity  w®*  of  the  c  frame  with  respect  to  the  i  fraune  is  simply  ncj. 

Thomson  (22:33-36)  derives  the  rotation  matrix  7l'‘  that  transforms  from  the  inertial 
frame  to  the  orbital  reference  frame.  The  process  of  rotating  from  the  i  frame  to  the  c 
frame  is  described  by  rotating  an  angle  ft  about  the  inertial  z  axis,  then  rotating  an  angle 
t  about  the  new  x  axis,  amd  finally  rotating  ui  angle  «  about  the  newest  z  axis.  The  amglee 
ft,  t,  «  are  often  referred  to  ats  a  3-1-3  Euler  angle  set.  The  rotation  matrix  is  formed 
by  the  product  of  three  intermediate  rotation  matrices: 

=  ^,(tt)7^,(t)7^,(ft) 
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where 


Tl^{v.)  = 


cos  u  sin  u  0 
—  sin  tt  cos  u  0 
0  0  1 


^x(t) 


10  0 
0  cos  t  sin  i 
0  —  sin  i  cos  i 


n,{ii)  = 


cos  n  sin  n  0 
—  sin  ft  cos  ft  0 
0  0  1 


The  resulting  rotation  matrix  is  calculated  to  be: 


cos  ft  cos  u  —  sin  ft  cos  t  sin  u 
—  cos  ft  sin  u  -  sin  ft  cos  *  cos  tt 
sin  ft  sin  t 


sin  ft  cos  tt  -f  cos  ft  cos  i  sin  tt  sin  i  sin  tt 

-  sin  ft  sin  tt  + cos  ft  cost  cos  tt  sintcostt 
-  cos  ft  sin  t  cos  t 


(2.2) 


Another  important  coordinate  transformation  is  that  relating  the  orbital  reference 
frame  to  the  body  frame.  The  rotation  matrix  we  derive  is  important  in  section  3.3 
where  an  expression  for  magnetically-induced  moment  is  derived.  Here  we  use  the  concept 
of  yaw,  roll,  and  pitch  angles  to  represent  the  rotational  angles  needed  to  rdate  the  frames. 
These  angles  are  collectively  referred  to  as  orientation  or  attitude  angles.  The  process  of 
rotating  from  the  c  frame  to  the  b  frame  is  described  by  rotating  am  angle  V'l  about  the 
£  axis,  then  rotating  an  angle  V*}  about  the  new  y  axis,  amd  Anally  rotating  an  angle  V's 
about  the  newest  z  axis.  Rirarott  (18:329-330)  refers  to  these  angles  with  this  particular 
sequaice  of  amgulau:  rotations  as  “Cardan  amides  of  the  Arst  kind”.  The  rotation  matrix 
TV*  is  formed  by  the  product  of  three  intermediate  rotaU:ion  matrices: 

=  Ti.{i>z)v^{^2)nMx) 
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where 


10  0 
0  cos  sin  V’l 

0  —  sin  cos 


'^»(V’2) 


cos  ^2  0  —  sin  V’a 
0  10 
sin  ^2  0  cos  ij)2 


cos  ^3  sin  03  0 
-sin 03  cos 03  0 
0  0  1 


The  resulting  rotation  matrix  is  calculated  to  be: 


cos  02  cos  0s  cos  0i  sin  03  +  sin  0i  sin  02  cos  03 
-  cos  02  sin  03  cos  0i  cos  08  -  sin  0i  sin  02  sin  03 
sin  02  -  sin  0i  cos  02 


sin  01  sin  03  -  cos  0i  sin  02  cos  03 
sin  01  cos  03  +  cos  0i  sin  02  sin  03 
cos  01  cos  02 


If  we  assume  that  the  orientation  angles  are  ‘'small”,  then  we  can  use  small-angle  assump¬ 
tions  to  linearize  the  71^*  matrix.  For  small  angle  0: 


sin^  »  0 

COS0  1 


The  products  of  small  angles  are  assumed  to  be  n^ligible  and  are  not  included  in  the 
linearized  form  of  the  matrix.  Thus  for  small  angles  0i,  02,  03,  the  linearized  form  of  1l‘* 


is: 


1 

-03 

02 


03 

1 

-01 


-02 

01 

1 


(2.3) 


We  assume  that  “small”  amgles  are  less  than  15*  for  the  purposes  of  this  research. 


Rotation  matrices  are  orthogonal.  The  inverse  of  an  orthogonal  matrix  is  equal  to 
its  tramspose.  Thus  it  is  easy  to  form  the  linearized  version  of  72*': 


72*'  =  72'=*'^  = 


1  -V’S  ^2 
^3  1  -A 


-A  A  1 


(2.4) 


2.2  Unforced  Equations  of  Motion 

The  tramslational  equations  of  motion  of  the  tracking  vehicle  relative  to  the  target 
vehicle  are  derived  by  Wiesel  (24:78-80)  and  by  Kaplan  (9:108-111).  The  derivation  of 
the  small-angle  relative  rotational  equations  of  motion  in  section  2.2.2  is  believed  to  be 
partially  original. 

2.2.J  Translational  Equations  The  complete  derivation  of  the  relative  translational 
equations  of  motion  is  not  shown  by  Wiesel  or  Kaplan.  All  steps  of  the  derivation  are 
explicitly  shown  here  to  illustrate  linearization  techniques  used  later  in  the  thesis. 

We  want  to  derive  an  expression  for  the  inertial  acceleration  of  the  tracking  vehicle 
and  expressions  for  the  forces  acting  on  it.  These  expressions  will  be  substituted  into 
Newton’s  second  law  to  give  the  translational  equations  of  motion: 

X;£  =  M,a  (2.5) 

where  is  the  sum  ail  external  forces,  Mt  is  the  mass  of  the  tracking  vehicle,  and  a  is 
the  inertial  accderation  of  the  tracking  vehicle. 

We  derive  the  inertial  acceleration  first.  Refer  to  the  position  vectors  shown  in 
figure  2.1.  The  inertial  position  vector  of  the  tradring  vehicle  is: 

£  =  lo  +  5  =  (»'o  +  x)ct  +  ycj  +  zci  (2.6) 
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where  is  the  distance  from  the  center  of  the  earth  to  the  c  frame  origin  and  y,  z  are 
the  c  frame  components  of  the  relative  position  vector  The  inertial  velocity  is: 


^  Ai 

=  Srr  +  w  xr 
at 

=  (ici  +  yc2  +  ics)  +  ncs  x  ((ro  +  x)ci  +  yc^  +  zcz) 

=  (x  -  ny)ci  +  (y  +  nro  +  na;)c2  +  ica 


The  inertial  acceleration  is: 


=  ((x  -  ny)ci  +  (y  +  nx)c2  +  zca)  +  ncj  x  ((x  -  ny)ci  +  (y  +  nro  +  nx)c2  +  zcs)) 

=  (x  -  2ny  -  n^fo  -  n^x)ci  +  (y  +  2nx  -  n^y)c2  +  zcs  (2.8) 


We  have  derived  the  accderation  part  of  Newton’s  second  law.  Now  •^e  decide  which 
external  forces  to  include  in  the  equations  of  motion.  There  are  several  types  of  external 
forces  that  could  act  on  the  tracking  vehide.  The  largest  of  these  in  magnitude  for  low 
earth  orbits  are  gravitational  force,  magnetic  force,  and  aerodynamic  drag.  In  section  3.3 
a  linearized  expression  for  ma^etic  force  is  derived.  Aerodynamic  drz^  is  described  in 
appendix  B.  Now  we  derive  a  linearized  expression  for  the  gravitational  force  The 
general  formula  for  the  gravitational  force  on  the  tracking  vehicle  is: 

Fj,  = - ^ — r 

— #  |.3 

Substituting  r  =  (ro  +  x)ci  +  yci  +  zcs  into  this  formula  yidds: 


_  -GM»M,((ro  +  x)ci  +  yc2  +  xcg) 


(rg  +  2roX  +  x*  +  y»  +  z*)l 
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The  tracking  vehicle  must  be  able  to  track  the  target  vehicle  at  “close”  range.  We 
assume  that  “close”  means  that  i  Tq,  j/  «C  ro,  and  z  <C  ro.  The  expression  for  gravita¬ 
tional  force  (eq.  2.9)  is  simplified  (linearized)  by  making  these  assumptions.  The  following 
steps  are  used  in  the  linearization  of  the  gravitational  force  expression: 

•  A  partial  binomial  series  expansion  is  used. 

•  Products  of  “small”  relative  position  components  r,  y,  and  z  are  neglected. 

The  binomial  formula  (19:110)  states  that  an  expression  of  the  form  (o  6)",  where  b  <  a, 
is  equivalent  to  the  series: 

(a  -b  6)"  =  a"  -1-  na^~^b  -!-•••  (higher  order  terms) 

Using  this  idea  to  simplify  the  denominator  of  equation  2.9: 

(r2-f2ror  +  x^-f  ~  (r^ -f  2rox)"^ 

~  To  *  -  3ro 

Substituting  this  result  into  equation  2.9  and  remembering  n  =  y/GM^Jr^ : 

Ff  ~  -GMQM,({ro  +  x)ci-\-yc2-\- zc3){ro^  -  rg*x) 

- r — (('■o  -  2x)ci  -I-  yc2  -I-  zca) 

^0 

~  -M,n^((ro  -  2x)ci  -f  yc^  -t-  zca)  (2.10) 
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Now,  substituting  the  inertial  acceleration  (eq.  2.8)  and  the  gravitational  force  (eq.  2.10) 
into  12  £  =  MtO,  we  obtain  the  relative  translational  equations  of  motion: 

x-2ny-Zn^x  =  (2-11) 

Mt 

y  +  2nx  =  (2.12) 

Mt 

z^-n^z  =  (2.13) 

Ml 

where  are  the  c  frame  components  of  the  sum  of  all  non-gravitational 

forces  acting  on  the  tracking  vehicle.  Wiesel  (24:78)  calls  this  set  of  equations  the  Clohessy- 
Wiltshire  equations.  These  differential  equations  ue  linear,  coupled,  and  have  constant 
coefficients.  When  £-other  —  we  call  these  equations  the  unforced  relative  translational 
equations  of  motion. 


2.2.1. 1  Solution  to  the  Unforced  Equations  The  unforced  equations  can  be 
solved  in  closed  form  (24:80-81)  to  yidd: 


x(0 

y(0 


z{t) 


/2  \  1  2 
-  (  “y(0)  +  3®(0)  j  cos  nt  +  "*(0)  sin  nt  +  4*(0)  +  — y(0) 

y(0)  -  (3y(0)  +  6n®(0))t  +  s*® 


-Im 

2(0)  cos  nt  +  —2(0)  sin  nt 
n 


(2.14) 

(2.15) 


(2.16) 


We  note  that  x{t)  contains  only  sinusoidal  and  constwt  terms.  z{t)  contains  purely  sinu- 
smdal  terms.  Interestingly,  y(t)  contains  the  linear  term  (3ir(0)  +  6nx(0))t  that  grows  with 
time  as  wdl  as  sinusoidal  terms  and  constant  terms.  This  means  that  the  tracking  vehicle 
will  move  steadily  away  from  the  target  vehide  in  the  y  direction  if  x(0)  and/or  y(0)  is 
non-zero.  We  see  in  section  2.2.2.1  how  the  soluticms  x(t),  y((),  z{t)  affiect  the  solutions  of 
the  rdative  small-angle  rotational  equations  of  motimi. 

The  unforced  rdative  translational  motion  of  the  tracking  vehicle  is  solely  dependent 
upon  its  initial  rdative  position  and  relative  vdodty.  Except  for  two  special  cases,  the 
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tracker  will  never  stay  at  its  initial  relative  position.  One  trivial  case  occurs  when  all  the 
relative  positions  and  velocities  are  zero.  The  origin  of  the  body  frame  and  the  origin  of 
the  orbital  reference  frame  are  then  coincident.  In  the  absence  of  any  perturbing  forces, 
the  tracker  will  stay  at  this  zero  reference  state.  A  more  interesting  special  case  occurs 
when  the  relative  velocities  are  all  zero,  the  relative  positions  a:(0),  z{0)  are  zero,  and  y(0) 
is  any  non-zero  value.  Then  the  tracker  will  remain  at  a  distance  j/(0)  ahead  or  behind  the 
target  vehicle  in  the  same  orbit. 


2. 2. 1.2  Steady-State  Behavior  Eventually  we  would  like  to  design  a  position 
controller  for  the  tracking  vehicle  which  will  allow  it  to  track  steady-state  values  of  a:,  y, 
and  z.  We  refer  to  these  steady-state  values  as  a;„,  y„,  and  z,,.  Now  we  examine  the 
steady-state  forces  required  to  maintain  steady-state  position  tracking.  In  the  steady-state, 
x  =  x  =  y  =  y  =  z  =  z  =  Q.  The  relative  translational  equations  of  motion  (eqs.  2.11, 
2.12,  2.13)  in  the  steady-state  become: 


Fy..  =  0 

Fi,,  —  Ti  MfZf, 


(2.17) 

(2.18) 
(2.19) 


So  the  steady-state  force  required  for  the  tracker  to  standoff  a  fixed  x  (vertical)  and/or 
z  (out-of-plane)  distance  from  the  target  vehicle  is  directly  proportional  to  the  standoff 
distance,  the  tracking  vehicle  mass,  and  the  square  of  the  orbital  mean  motion.  We  note 
that  three  times  as  much  force  is  required  to  maintain  the  x  component  of  the  standoff 
distance  as  is  required  to  maintain  the  z  component.  This  fact  makes  out-of-plane  standoff 
less  ‘^costly”  than  vertical  in-plane  standoff.  No  steady-state  force  is  required  to  maintmn 
the  y  component  of  the  standoff  distance.  We  note  that  a  force  in  the  negative  x  direction 
(down)  is  needed  to  maintain  a  constant  positive  (up)  x  component  of  the  standoff  distance. 
We  also  observe  that  the  required  steady-state  force  for  a  given  standoff  distance  becomes 
smaller  as  the  reference  orbital  radius  ro  increases  since  =  GM^/rg. 


2.2.2  Rotational  Equations  Kaplan  (9:201-202)  and  Wiesel  (24:145-148)  derive  the 
small-angle  rotational  equations  of  motion  of  the  gravity-gradient  satellite  for  the  purpose 
of  studying  attitude  stability.  The  body  attitude  angles  that  these  authors  use  are  defined 
with  respect  to  an  orbital  reference  frame  attached  to  the  body.  The  important  difference 
in  the  following  derivation  is  that  the  orientation  angles  are  defined  with  respect  to  an 
orbital  reference  frame  attached  to  another  body  (the  target  vehicle).  The  equations  we 
derive  are  the  rotational  equivalent  of  the  Clohessy- Wiltshire  translational  equations  of 
motion. 

We  start  with  Euler’s  general  rotational  equations  of  motion: 

AtJi  -|-  (C"  —  ~  hfi 

Boi^  "I"  (■‘^  —  C)witj3  =  M2  (2.20) 

CcJa  -|-  (5  —  A)(jJi(jJ2  =  M3 

where  A,  B,  C  are  the  principal  moments  of  inertia,  uji,  u>2,  ^3  are  the  body  frame 
components  of  the  inertial  angular  velocity  vector  <£,  and  Mi,  M2 ,  M3  are  the  body  frame 
components  of  the  sum  of  all  external  moments.  We  need  to  express  the  inertial  angular 
velocity  u  and  the  external  moment  M  in  terms  of  the  orientation  angles  V'l)  V'2)  and  V's- 

There  are  several  possible  external  moment  sources  acting  on  the  tracking  vehicle. 
Only  the  gravitational  moment  M^  is  examined  in  this  section.  In  section  3.3  an  expression 
for  magnetic  moment  is  derived.  Any  moment  au’ising  from  aerodynamic  drag  is  assumed 
to  be  negligible. 

Wiesel  (24:145-147)  derives  the  gravitational  moment  in  body  framie  coordinates: 


M,.  = 

^^YZ{C  -  B) 

= 

^^XZiA  -  C) 

(2.21) 

II 

—p^XY{B  -  .4) 

where  X,Y,  Z  are  the  body  frame  components  of  the  position  vector  r.  Here  the  derivation 
diverges  from  that  of  Wiesel  and  Kaplan  because  our  position  vector  r  includes  the  relative 
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position  components  x,  y,  z.  We  need  to  express  the  given  gravitational  moment  in  terms 
of  the  orientation  angles.  We  begin  by  expressing  r  in  body  frame  coordinates: 

*r  = 

where  the  linearized  form  of  72®*  is  used  (eq.  2.3)  because  the  orientation  angles  are  assumed 
to  be  small.  Neglecting  the  products  of  small  variables  x,  y,  z,  ipi,  ^’25  a^nd  V'a  we  get: 

^  (ro  +  a:)ii  +  (-V»3»‘o  +  y)b2  +  (V’2»‘o  +  2)^3 

which  can  also  be  stated: 

*r  =  Xbi  +  Yb2  +  Zbs 

where 


X  a  To  +  ® 
y  ~  -^3To  +  y 
Z  ~  ^2»*0  +  2 

Again  ne^ecting  the  products  of  small  variables,  the  terms  YZ,  XZ,  XY  are  evaluated: 

YZ  ~  0 

XZ  ~  rlxIfi  +  Toz  (2.22) 

XY  ~  -rli>3  +  Toy 


The  denominator  r~*  of  the  moment  terms  is  linearized: 

r~®  =  (fo  +  2ro*  +  **  +  y*  +  2*)“  * 

~  (fo  +  2roa:)"*  (2.23) 

~  To  ®  -  5ro  ®ar 
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Finally  substituting  equations  2.22  and  2.23  into  the  gravitational  moment  equations 
(eq.  2.21),  we  get  the  linearized  form  of  M^: 


~  0 

~  3n"(V»2  +  -)(A-C) 

^0 

M,,  2;  3n*(-^3  +  ^){B  -  A) 
To 


We  now  have  the  gravitational  moment  part  of  the  small-angle  rotational  equations 
of  motion.  Now  we  derive  the  inertial  ang^ilar  velocity: 


u  =  nca  -1-  V’lhi  -|-  +  V»3^3 

,  ■  C3  -1-  ^ibi  -h  ^363 


where  the  linearized  form  of  y/GM^Jr^  is: 


=  (GM,)*(r»  +  2r,»+.’  +  y’+z>)-! 

Then  the  inertial  angular  velocity  becomes: 

u>  =  (GMe)i(ro  *  -  **)cs  -I-  ^iSi  -I-  v>3ftj  +  ^3*3 

Using  the  linearized  rotaticm  matrix  H"*  (eq.  2.3)  to  transform  the  £3  term  into  the  S 
frame  (again  n^ecting  the  products  of  small  variables)  we  derive  the  linearized  form  of 
the  inertial  angular  velocity: 


•  A  •  A  3fix 

w  =  (-71^2  +  ^i)bi  +  (»V»1  +  ^)b2  +  (n-h^3-  —)b3 


(2.24) 
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We  use  equation  2.24  to  find  oii,  tJ2,  ^3: 


^  .A  ••  •  ^  ••  ^713^  A 

'-t:^  =  {i>i  -  nil)2)bi  +  (t^2  +  nipi)b2  +  (^3  -  — )^>3 
ut  ZTq 

and  the  products  a;2W3,  ujiu>2: 

*  *  3yix  * 

W2W3  =  {ip2  +  nij}i)(ip3  +  n  -  - — )  ~  TiV’2  +  nVi 

27*0 

*  ’  ^TtCC  * 

(Jia;3  =  (V^i -«V’2)(03  +  n- — )  7101-0^^2 

zro 

^1072  =  (V’l  -  raV’2)(^2  +  nV»l)  ~  0 

Substituting  these  results  into  the  general  rotational  equations  of  motion  (eq.  2.20),  we 
obtain  the  small-angle  relative  rotational  equations  of  motion: 

^(^1  -  77^2)  +  {C  -  B){ni)2  +  =  0  (2.25) 

B(^2  +  n-i/ji)  -f  (A  -  C)(n^i  -  471^^2 - )  =  0  (2.26) 

To 

C(^3-^)  +  {B-A)i3n^tP3-^)  =  0  (2.27) 

These  differential  equations  are  linear,  constant-coefficient,  and  coupled.  The  -3n^z/ro 
term  in  the  second  equation  and  the  —3nx/2ro  and  —3n^ylro  terms  in  the  third  equation 
do  not  appear  in  the  rotational  equations  of  motions  derived  by  Kaplan  and  Wiesel.  These 
“extra”  terms  arise  from  the  definition  of  the  body  orientation  angles.  These  angles  were 
defined  with  respect  to  an  orbital  reference  frame  attached  to  miother  body  (the  target 
vehicle)  instead  of  an  orbital  reference  frame  attached  to  the  tracking  vehicle. 


2.2.2. 1  Unforced  Equation  Observations  The  solutions  x(t),  y{t),  z{t)  (eqs. 
2.14,  2.15,  2.16)  to  the  relative  translational  equations  of  motion  can  be  inserted  into  the 
“extra”  terms  in  the  relative  rotational  equations  of  motion  to  yield: 


A(^i  -  71^2)  +  (C  -  B){nrj)2  -|-  =  0 

;  ,,  37»^(z(0)  COS  Ttf ^i(O)  sin  7»f)  v 

+  nijji)  -f  (^  -  C){nrl>i  -  47i’V'2 - ^ - — - -)  =  0 


»"o 
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C(fe  _  Mi2»(0)  ±3n^m  sinnl  +  i(Q)  cosjl)  ^ 

2ro 

3n^(y(0)  -  (3y(0)  +  6na;(0))f  +  (^y(O)  +  6a:(0))  sin  nt  +  f  ^(0)  cos  nt  -  f  i(0))  _  ^ 

To 

Note  that  these  equations  become  the  equations  shown  in  Wiesel’s  text  (24:148)  when  the 
relative  positions  and  velocities  are  zero.  We  call  the  equations  that  we  have  derived  here 
the  “unforced”  equations  because  we  assume  no  forces  other  that  gravitational  are  acting 
on  the  tracking  vehicle.  In  another  sense  we  could  think  of  these  equations  as  the  rotational 
equations  shown  in  Wiesel’s  text  with  added  forcing  functions  which  result  from  non-zero 
initial  relative  positions  and  velocities.  Wiesel  shows  that  the  solutions  to  the  unforced 
equations  are  oscillatoiy  in  nature  as  long  as  the  condition  C  >  B  >  A  is  satisfied. 

Although  we  do  not  explicitly  solve  the  relative  rotational  equations  of  motion  for 
V’iCOj  V’3(0>  make  general  comments  about  the  nature  of  the  solutions. 

We  assume  that  C  >  B  >  A.  Then  ^i(t)  is  purely  sinusoidal.  V’aCO  is  ^dso  sinusoidal 
because  the  forcing  function  z(t)  is  sinusoidal.  V’sCO  i®  more  interesting  because  while 
x{t)  is  sinusoidal,  y(t)  contains  linear  and  constant  terms  in  addition  to  sinusoidal  terms. 
Therefore  if  a:(0)  is  non-zero,  V’3(0  grows  larger  with  time  until  the  small-angle  assumption 
becomes  invalid.  Then  the  linearized  equation  in  V’aCO  ^^ot  apply  to  the  non-linear 
behavior  of  V’aCO- 

2.2.S.S  Steady  State  Behavior  We  want  to  eventually  design  an  attitude  con¬ 
trol  system  for  the  tracking  vehicle  which  will  allow  it  to  regulate  the  steady-state  values 
of  V'l,  ^27  and  ^3  to  zero.  We  refer  to  these  steady-state  values  2is  ^2,, 7  and  V’3,,.  Now 
we  examine  what  moments  are  required  to  maintain  steady-state  attitude  regulation.  In 
the  steady-state,  all  the  orientation  angle  rates  and  position  rates  go  to  zero.  Note  that  the 
relative  position  variables  are  allowed  to  be  non-zero  in  the  steady  state.  The  small-angle 
relative  rotational  equations  of  motion  (eqs.  2.25,  2.26,  2.27)  in  the  steady-state  become: 

Mr..  =  n\C-B)t/;,„  (2.28) 

My..  =  -4n^(A-C)V>2..-^"'^--  (2.29) 

To 
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M, 


(2.30) 


=  3n^(B  -  A)‘i}}z„  - 


Zn\B  -  A) 


^0 


Vss 


We  usually  want  V’l,,)  V’2„»  stud  rj)z„  to  be  equal  to  zero  to  avoid  rotational  inertia  during 
a  docking  maneuver.  From  these  equations  we  then  expect  that  Mr,,  is  equal  to  zero, 
is  constant  and  proportional  to  the  steady-state  z  component  of  position,  and  is 
constant  and  proportional  to  the  steady-state  y  component  of  position.  We  expect  that  the 
forces  required  to  generate  these  steady-state  moments  are  much  smaller  than  the  forces 
required  for  steaxiy-state  position  standoff  when  and  V’a,,  are  zero. 
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III.  Magnetic  Forces  and  Moments 


In  this  chapter  we  complete  the  linearized  relative  translational  and  rotational  equa¬ 
tions  of  motion  by  deriving  expressions  for  magnetic  force  and  moment. 


3. 1  Geomagnetic  Field  Description 

The  geomagnetic  field  strength  B  is  equal  to  the  negative  gradient  of  the  scalar 
magnetic  potential  (14:2.22): 

5  =  -  V 

The  #m  that  one  chooses  for  analysis  purposes  depends  upon  the  level  of  accuracy  one 
wants  and  the  complexity  that  one  is  willing  to  tolerate.  Unfortunately^ the  more  accurate 
of  a  model  desired,  the  more  complex  becomes.  The  most  accurate  model  is  the 
multipole  expansion  model  (also  known  as  the  spherical  harmonic  analysis  model)  (7:41). 
The  scalar  magnetic  potential  for  this  model  is: 

=  r®  X)(— )’*'^^  Yli9n  cosmA  -I-  h”  sin mX)P;^ (sin 6)  (3.1) 

B=1  ”  m=0 

where  r®  is  the  radius  of  the  earth,  6  is  the  magnetic  co-latitude,  A  is  longitude,  smd 
/^(sini)  is  the  associated  Legendre  function  of  order  m  and  degree  n  (19:149).  The 
coefficients  g^  and  are  known  as  the  Schmidt  coefficients.  Several  are  tabulated  in 
NASA  TM  82478  (14:2.23). 

We  will  instead  use  the  simple  dipole  model  (21:33-34)  to  represent  the  geomagnetic 
fi^d.  The  scalar  magnetic  potential  for  this  model  is  equivalent  to  the  first  term  of  the 
scalv  magnetic  potential  expression  of  the  multipcde  expansion  model  (eq.  3.1).  The 
farther  one  moves  from  the  earth,  the  more  accurate  the  simple  dipole  model  becomes 
because  the  higher  order  multipole  terms  decrease  in  magnitude  more  quickly  tham  the 
dominant  dipole  term.  The  scalar  magnetic  potential  for  the  simple  dipole  model  is: 


$ 


m 


M'L 

r* 


(3.2) 
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X  =  lll 


g3 


Figure  3.1.  Magnetic  Dipole  Moment  Vector  in  the  Greenwich  Frame 

where  m  denotes  the  magnetic  dipole  moment  of  the  earth  and  r  is  the  position  vector 
from  the  center  of  the  earth  to  the  tracking  vehicle.  The  m  vector  originates  at  the  center 
of  the  earth  and  passes  through  the  earth’s  surface  at  the  austral  magnetic  pole  in  the 
southern  hemisphere.  This  modd  may  give  a  magnetic  held  magnitude  error  of  up  to  30% 
at  the  earth’s  surface  according  to  Hess  (7:39). 

S.2  Transformation  to  the  OrhittA  Reference  Frame 

The  rdative  translaticMial  equations  of  motion  (eqs.  2.11,  2.12,  2.13)  were  derived  in 
orbitjd  reference  frame  (c)  coordinates.  An  expression  for  the  magnetic  force  ^  in  c  frame 
coordinates  is  derived  in  this  section.  This  allows  the  convenient  derivation  of  the  magnetic 
forces  acting  on  the  tracking  vehicle  through  the  use  of  the  Lorentz  equation  (eq.  1.1). 
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Using  the  simple  dipole  model,  the  geomagnetic  field  is  expressed: 


B  =  -V#, 


In  order  to  derive  we  must: 

•  Express  m  in  the  c  frame  (^m)- 

•  Perform  the  dot  product  m  •  r. 

•  Apply  the  gradient  operator  V- 

This  particular  method  of  deriving  the  magnetic  field  in  terms  of  orbital  reference  frame 
coordinates  was  suggested  by  Wiesel  (25). 

The  magnetic  dipole  moment  vector  is  easily  visualized  in  a  coordinate  frame  we 
designate  the  greenwich  frame  (see  figure  3.1).  The  g  frame  differs  from  the  inertial  geo¬ 
centric  equatorial  frame  (i)  by  a  single  rotation  about  the  inertial  z-axis  (ta).  The  angle 
of  rotation  Og  locates  the  greenwich  prime  meridian  with  respect  to  the  inertial  x-axis  (tj). 
This  angle  is  called  the  sidereal  time  and  is  calculated  using  the  following  simple  formula 
(1:99): 

9,  =  9g,+uj9it-to)  (3.4) 

where  9,0  is  the  value  of  at  reference  time  to  a^d  is  the  angular  vdocity  of  the  earth’s 
rotation.  The  American  Ephemeris  and  Nautical  Almamac  can  be  used  to  determine  . 
The  angles  S  and  A  shown  in  figure  3.1  ve  the  latitude  and  longitude,  respectively,  of  the 
austral  magnetic  pole.  These  an^es  completely  determine  the  direction  of  the  magnetic 
dipcde  moment  vector  si  in  the  greenwich  frame: 

=  m  cos  S  cos  Xgi  -I-  m  cos  S  sin  -I-  m  sin  Sgg  (3.5) 

where  the  magnitude  of  the  magnetic  dipole  moment  m  is  aq>proximatd[y  8.1  •  10^'’  tesla- 
meter®.  For  the  austral  pole  6  ~  —78.5*  and  A  ~  111*. 
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The  transformation  of  m  from  the  g  frame  to  the  c  frame  yields  the  c  frame  components 
of  m; 

mi  =  m(cos  6  cos  A(cos  (fi  —  cos  «  —  sin  (ft  -  )  cos  i  sin  u) 

+  cos  6  sin  A(8in  (ft  —  cos  u  +  cos  (ft  —  )  cos  i  sin  u) 

+  sin  6  sin  i  sin  u) 

m2  =  m(cos^  cos  A(— cos(ft  —  ^,)sin«  —  8in(ft  -  ^>,)cos  t  cosu) 

+  cos  6  sin  A(-  sin  (ft  —  sin  it  +  cos  (ft  —  (?,)  cos  i  cos  u) 

+  sin  6  sin  i  cos  it) 

m3  =  m(cos  6  cos  A  sin  (ft  —  9^)  sin  i 
-  cos  S  sin  A  cos  (ft  —  d,)  sin  t 
+  sin  S  cos  i) 

which  when  simplified  by  the  trigonometric  identities: 

sin(i4±5)  =  sin  j4  cos  ±  cos  A  sin  5 

cos  (j4  ±  5)  =  cos  A  cos  ^  ^  sin  A  sin  B  (3.6) 

become: 

mi  =  m(co8^  cositco8(ft  -  —  A)  —  co8tco8^sinu8in(ft  -  -  A)  +  sintsintfsintt) 

m2  =  m(— cos  5  sin  It  cos  (ft  —  —  A)  —  cos*  cos  ^  cos  It  sin  (ft  —  -  A)  +  sin  t  sin  ^  cos  It) 

m3  =  m(sin  I  cos  ^8ia(ft  —  —  A)  +  cos  *  sin  #) 
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We  introduce  some  simplifying  notation  before  performing  the  dot  product  (m-r). 
Let 

a  =  cos^cosucos(fl  —  —  A)  —  cos*  cos^sinusin(f2  —  —  A)  +  sintsinrfsinu 

/?  =  —  cos  S  sin  «  cos  (Q  —  —  X)  —  cos  *  cos  S  cos  u  sin  (fl  —  —  A)  +  sin  *  sin  S  cos  u 

7  =  sin*cos^sin(ft  —  —  A)  +  cosisin^ 

where  we  note  for  future  reference  that  a,  (3  are  composed  of  purdy  sinusoidal  components 
while  7  indudes  the  constant  term  cos  *  sin^  as  well  as  a  sinusoidal  component. 

Now  we  can  write 'm  =  m(aci  +  /3c2  +  7C3).  Remembering: 

‘r  =  (ro  +  x)ci  +  yc2  +  ZC3 


the  dot  product  m  •  r  becomes: 

rn-r  =  m(o(ro  +  *)  +  +  7^) 


The  remaining  hurdle  in  the  evaluation  of  the  magnetic  field  B  in  c  frame  coordinates 
is  the  application  of  the  gradient  operator: 


_  (  (ai 


(rp  +  j)  +  +  iz) 


2roX  +  **  +  p*  +  x' 


i)i) 


which  can  also  be  expressed: 


D_  f  ^  ^  ^  \  (  (ar(ro  +  x)  +  /3y  +  jz)  \ 


Performing  the  partisd  differentiation,  making  a  binomial  approximation  to  the  resulting 
denominators,  and  neglecting  the  products  of  the  small  relative  position  components  x,  p, 
z,  we  m*rive  at  the  linearized  expression  for  the  magnetic  field  in  orbital  reference  frame 
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bl  (up) 


Figure  3.2.  Design  #  1  Conductor  Configuration 


coordinates: 

=  ^((2aro  -  6aa:  +  3/9y  +  37z)ci  +  ( -Btq  +  3/3®  +  3tty)c2  +  (-7ro  +  37®  +  3az)c3)  (3.7) 

Appendix  A  shows  plots  of  these  orbital  frame  components  of  the  magnetic  field  as  func¬ 
tions  of  time  assuming  the  relative  position  components  *,  p,  z  are  all  zero.  The  plots 
are  dtme  for  a  400  km  orbit  at  several  different  inclination  angles.  From  these  plots  the 
reader  f  n  see  that  components  disappeau'  at  certain  times  in  the  orbit.  Both  in-plane 
components  or  the  out-of-plane  component  disappear  watirely  at  some  inclinations.  This 
becomes  important  in  the  next  chapter  where  we  investigate  the  subject  of  controllability. 

S.S  Magnetic  Force*  and  Moments  for  Design  #  1 

An  dectrical  conductor  configuration  is  desired  that  allows  independent  control  over 
every  component  of  tracking  vehicle  position  and  attitude.  Two  designs  are  proposed  in 
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this  section  to  accomplish  that  goal.  In  chapter  four  we  show  that  the  second  design  is 
more  controllable  than  the  first  design.  Design  #  1  is  shown  in  figure  3.2.  There  are  a  total 
of  six  conductors:  two  on  the  in-plane  body  x-axis  (one  “up”  and  one  “down”),  two  on  the 
in-plane  body  y-axis  (one  “east”  and  one  “west”),  and  two  on  the  out-of- plane  z-axis  (one 
“north”  and  one  “south”).  Lightweight  structural  webbing  exists  between  the  conductors 
to  provide  rigidity.  The  fiUed-in  circles  at  the  ends  of  the  conductors  shown  in  figure  3.2 
represent  the  plasma  generators  (PMGs)  mentioned  in  chapter  one. 

S.3.1  Magnetic  Forces  The  conductor  lengths  and  currents  are  labeled  in  figure 
3.2.  Each  conductor  current  is  controlled  independently  of  the  other  currents.  The  Lorentz 
equation  (eq.  1.1)  is  used  to  determine  the  force  on  each  conductor: 


Fg  =  ILxB  (3.8) 

We  remind  ourselves  that  the  length  vector  L  is  defined  to  point  in  the  direction  of  positive 
current  flow.  The  individual  forces  on  each  conductor  are  summed  to  determine  the  overall 
magnetic  force  on  the  tracking  vehicle: 


The  summation  IL  is  conveniently  expressed  in  the  body  frame: 


^2^*—  ~  +  ^3.13^3  +  Izjahs 


—  +  ^4)^1  +  hih,  +  +  h.)b3 
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Using  the  linearized  form  of  the  rotation  matrix  72*®  (eq.  2.4)  to  transform  to  the  c  frame 
we  get: 


+ -fli)  ~  4-  V’2^3(f3.  +  f3.))ci 

+  (i’shiJi,  +  hi)  +  hih,  +  /2»)  —  i’lUih,  +  h,))c2 
+  (— ■02^i(^i.  +  Ju)  +  +  hj)  +  IsC^s.  +  h,))c3 


We  derived  an  expression  for  in  the  previous  section  where  we  linearized  with 
respect  to  the  “small”  position  and  attitude  variables  x,  y,  2,  V’2>  and  ^3.  In  the  next 

chapter  we  use  control  terminology  to  define  these  variables  as  the  states  of  the  tracking 
vehicle  system.  When  we  perform  the  cross  product  in  the  present  derivation  we  will 
linearize  with  respect  to  the  controls  as  well  as  with  respect  to  the  states.  The  tracking 
vehicle  controls  are  the  conductor  currents.  The  products  of  states  and  the  products  of 
states  and  controls  are  neglected.  This  assumption  greatly  simplifies  the  cross  product 
^£s  =  X  We  calculate  the  total  linearized  magnetic  force  on  Design  #  1: 


ip 


=  m 


\ 


4- 


”0 


Cl 


4-  m 


4-  m 


^20/3(73,  4- 13.)  7^i(-fi.  4-  hi) 


r3 

”0 


)c. 


(3.9) 


(- 


l^h{h,  4-  hi)  20/2(72.  4-  l2„y 


rg 


C3 


3.S.2  Magnetic  Moments  We  use  the  following  general  expression  to  evaluate  the 
magnetic  moment  on  the  tracking  vehicle: 


where  *r  r^resents  the  moment  arm  from  the  tracking  vehicle  center  of  mass  to  the  center 
of  applied  ma^etic  force  for  each  conductor.  We  assume  that  the  conductors  are  short 
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enough  so  that  we  can  reasonably  assume  that  the  center  of  magnetic  force  of  a  conductor 
coincides  with  the  geometric  center  of  that  conductor. 

The  magnetic  field  in  c  frame  coordinates  is  known  (eq.  3.7).  We  want  to  calculate 
the  magnetic  moment  in  terms  of  body  frame  coordinates  in  order  to  easily  insert  the 
resulting  expression  into  the  small-angle  relative  rotational  equations  of  motion.  Using  the 
linearized  form  of  71®*  (eq.  2.3)  to  transform  B_  to  the  b  frame  yields: 

i  m  « 

*5  =  -r((2Q:ro  -ax  +  ZjSy  -f-  372  -1-  7t/»2  -  /Jt/»3)6i 

^•a 

+  (-0ra  +  3/3x  -f-  3ay  -  'yrpi  -  2q;03)62  (3.10) 

+  (-7^0  +  37*  +  3az  +  /?V’i  +  2atj}2)bz) 

where  again  we  have  ne^ected  the  products  of  relative  position  and  attitude  components. 
Then  *M  is  expressed: 

X  (A./iSi  X  *B)  -  X  (/i,/i6i  X  *5) 

+  I&2  X  {hMk  X  *5)  -  |S2  X  {hjzk  X  '’B) 

+  1^3  X  {hMz  X  *5)  -  1^3  X  (/3./363  X  *£) 

which  simplifies  to: 

‘M  =  |6ix((/i.-/,>ix‘£) 

+  ^^2  X  ((/2,  —  Izjibi  X  *£) 

+  1*3  X  ((/a.  -  h,%  X  *5) 
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Performing  the  cross  products  and  linearizing  with  respect  to  the  states  and  the  controls 
we  find  the  total  linearized  magnetic  moment  acting  on  the  tracking  vehicle: 


+  (^(A.  -  A.) +  -/,))«, 


(3.11) 


3.3.3  Complete  Equations  of  Motion  Now  that  we  possess  expressions  for  the  total 
magnetic  force  (eq.  3.9)  and  moment  (eq.  3.11),  we  can  write  the  complete  relative  trans¬ 
lational  and  rotational  equations  of  motion  for  Design  #  1.  The  translational  equations  of 
motion  are: 


X  —  2ny  —  3n^x 
y  -|-  2nx 
z  -|-  n^z 


m'jliih.  +  hj  +  h.) 

Mtrl  M,tI 

2mal3{h^  +  h.)  rn-ylijh,  +  Ii,) 
Mtr%  Mtrl 

+  A^)  2ma/2(/2.  + /2„) 

M.rg  MiTf 


(3.12) 

(3.13) 

(3.14) 


The  rotational  equations  of  motion  are: 


;•  n(-A  -B  +  C)  ;  n'^(C  -  B)  , 

A  +  - /  -  -  — hi  = 


mall  r  N  j  s 


Arl 


Arl 


B 


B  ,  5ro 

'^(T  -T  -T  ) 


••  3n  .  37*2(5  -  A)  ,  37*2(5  -  A) 


(3.15) 


(3.16) 


(3.17) 


where  a,  /?,  and  7  were  defined  on  page  3-6.  We  observe  that  these  equations  are  highly 
coupled  with  respect  to  the  controls  as  well  as  the  states.  We  also  note  that  the  and  ^^>2 
equations  for  this  design  have  purely  sinusoidal  forcing  functions  on  the  right  hand  side. 
When  these  terms  become  zero  in  either  of  these  equations,  “control  authority”  has  been 
lost  in  that  equation.  This  behavior  is  a  first  hint  that  Design#  1  may  be  inferior. 


3.4  Magnetic  Forces  and  Moments  for  Design  if  Z 

Design  #  2  is  shown  in  figure  3.3.  There  are  a  total  of  eight  conductors:  two  01 
the  in-plane  body  x-axis,  two  on  the  in-plane  body  y-axis,  two  that  form  the  out-of-plane 
northern  cross,  and  another  two  that  form  the  southern  cross. 
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3.4.1  Magnetic  Forces  Each  conductor  length  and  current  is  labeled  in  figure  3.3. 
The  Lorentz  equation  (eq.  1.1)  is  used  again  to  determine  the  force  cm  eadi  conductor: 

Fb=  ILxB  (3.18) 

The  individual  forces  on  each  conductor  are  summed  to  determine  the  overall  magnetic 
force  on  the  tracking  vehicle: 

The  summation  IL  is  conveniently  expressed  in  the  body  frame: 

+  /2./2^2  + 

=  (^(A,  +  ht)  +  ^ei(A.  +  A.))ftl  +  {hih,  +  hy.)  +  le,(^2.  +  h.))b2 
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Using  the  linearized  form  of  the  rotation  matrix  7^**  (eq.  2.4)  to  transform  the  current- 
length  product  to  the  c  frame  we  get: 


=  (/l(/u+/lJ  +  /c.(/u  +  /l.)-^3/2(/2.  +  /2j-V’3U/2.+/2.))^l 

+  (V’3^(^l.  +  lu)  +  +  A.)  +  hih.  +  ^2.)  +  +  f2,))c2 

+  (  — V’2^l(-fl.  +  /l^)  —  V’2^Ci(f^l»  +  A.)  +  +  ■f2„)  +  +  ^2.))C3 


Performing  the  cross  product  X  *:fi)  a«d  linearizing  we  calculate  the  total 

linearized  magnetic  force  on  Design  #  2: 


=  m 


'rhjh.  +  /2.) 


(  ^3  J.3  J  *^1 

^/7/i(/i!  +  /u)  ,  7/c.(/u  +  /i.)V, 

V  To 

m 


(3.19) 


fihjhn  +  ^i<)  ^hiih.  4-  h.)  2al2{l2,  +  h^)  2alcj{l2.  +  h.T 


C3 


3.4.2  Magnetic  Moments  As  for  Design  1  we  use  the  following  general  expression 
to  evaluate  the  magnetic  moment  on  the  tracking  vehicle: 

=  Yj(*r  X  (/*X  X  *B)) 

where  *B_  was  calculated  in  equation  3.10.  Then  ^M_  is  expressed: 

=  |iiX(/,./Ax‘£)-|6,  x(AMx*5) 

+  1^2  X  (/,./,&,  X  *£)  -  X  {I2j2b2  X  ‘£) 

+  lih  X  {hJcM  X  *£)  +  /,fr3  X  {I2JJ2  X  ‘B) 

-  lih  X  (A.lc.Si  X  *£)  -  1,63  X  X  »B) 
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which  simplifies  to: 


‘M  =  X  ((/i.  - /i  Jfei  X  *5) 

+  1*2  X  ((/,.  -  hjk  X  '*B) 
+  lilc^k  X  ((/i.  -  li.)k  X  *5) 

+  /.U3X((J2.-/2.)62X»5) 


Performing  the  cross  products  and  linearizing  with  respect  to  the  states  and  the  controls 
we  find  the  total  linearized  magnetic  moment  acting  on  the  tracking  vehicle: 


‘M  = 


(3.20) 


3.4.3  Complete  Equations  of  Motion  Now  that  we  possess  expressions  for  the  total 
magnetic  force  (eq.  3.9)  and  moment  (eq.  3.11),  we  can  write  the  complete  relative  trans¬ 
lational  and  rotational  equations  of  motion  for  Design  #  2.  The  translational  equations  of 
motion  are: 


X  -  2ny  -  3n’x 


y  -I-  2ni 
z  +  n*2 


myltih,  +  /j.) 

+  h.) 

M,rg 

myliili^  +  lu) 

"»7^.(A.  +  /i.) 

M,r? 

M.rg 

4-  /14) 

m^/,.(/u  +  /i.) 

Af.r* 

M,rg 

2ma/j(/,.  -1-  /j. 

,)  2ma/c,(/2, -I-/2.) 

M,r? 

MfTQ 

The  rotational  equations  of  motion  are: 


(3.21) 

(3.22) 

(3.23) 


(3.24) 
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mall 


""“•2/r  T  \  f  J  r  \ 


..  njA  +  B-C)  _  4n^A-C)  _  3n  V  -  C)  ^ 

^  B  5  ^ro 

m/3/Lr  r  ^ 


2Br3 


(A.-AJ- 


5r3 


(/2.-I2.) 


;•  3n  .  3n\B  -  A)  ,  3n\B  -  A) 

— k-^'= 


(3.25) 


(3.26) 


where  a,  /?,  and  7  were  defined  on  page  3-6.  We  note  the  presence  of  7  terms  in  the  ipi 
and  ^2  equations.  Unlike  Design^  1,  control  authority  in  low  inclination  orbits  does  not 
appear  to  be  lost  in  these  equations  due  to  the  constant  term  contained  in  the  definition 
of  7. 
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IV.  Tracking  and  Control  Issues 


This  chapter  explores  the  issue  of  controller  design  for  the  tracking  vehicle.  The 
controller  is  designed  to  allow  the  tracking  vehicle  to  stay  a  fixed  distance  away  from  and 
maintain  a  constant  attitude  with  respect  to  a  target  vehicle.  Design  #  2  is  shown  to  be 
more  controllable  at  certain  points  in  low-inclination  orbits  tham  Design  #1.  A  MATLAB 
(16)  program  (CONTROLLER.M  -  see  appendix  F)  is  written  for  controller  design  and 
evaluation.  Also  a  new  control  theory  concept  we  call  steady-state  controllability  is  defined 
and  explored. 

4.1  State  Space  Representation  of  Equations  of  Motion 

We  use  a  state  space  representation  to  model  the  POTV’s  electrodynamic  interaction 
with  the  geomagnetic  field.  The  state  space  model  lends  itself  to  the  understamding  of 
system  behavior  and  controller  design.  Our  linearized  system  is  time- varying.  The  state 
space  model  takes  the  form  (3:363): 

i  =  Fx  -1-  G{t)u 
y  =  Hx 

where  the  state  vector  x  consists  of  the  relative  positions,  attitude  angles,  and  their  rates. 
The  number  of  states  in  a  system  is  the  system  order  N.  Our  system  is  twelfth-order. 

*  =  X  X  y  y  z  z  tl>2  i’l  '/’3  i’s 

The  output  vector  y  consists  of  the  system  variables  we  are  interested  in  controlling,  namely, 
the  relative  position  and  attitude  states.  The  dimension  of  y  is  6  x  1. 

y  =  ^  a:  y  z  ^2  V’s 
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The  control  vector  u  consists  of  the  conductor  currents.  The  dimension  of  for  Design  # 
1  is  6  X  1.  The  dimension  of  ttj  for  Design  #  2  is  8  x  1. 


F  is  the  12  x  12  system  matrix.  This  matrix  is  the  same  for  both  designs.  The  solution 
x{t)  of  i  =  Fx  represents  the  unforced  behavior  of  the  system.  G{t)  is  the  time-varying 
input  matrix.  It  shows  how  the  controls  enter  the  system  to  affect  the  states.  For  Design 
#  1,  <j  is  12  X  6.  For  Design  #  2,  (7  is  12  x  8.  H  is  the  6  x  12  output  matrix.  H  selects 
the  states  we  wish  to  control. 

The  six  second  order  differential  equations  of  motion  for  each  design  (eqs.  3.12-3.17 
and  eqs.  3.21-3.26)  must  be  rewritten  as  a  system  of  twelve  first  order  equations  in  order  to 
use  the  state  space  model.  This  process  is  described  by  Waltman  (23:2,88).  The  resulting 


state  space  matrices  are: 

the  system  matrix  F: 
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i(t)  for  Design  #  1: 
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the  input  matrix  GjCO  for  Design  #  2: 
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and  the  output  matrix  H-. 
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Note  that  the  rank  of  the  input  matrices  cam  never  be  greater  than  five. 

Six  of  the  rows  aure  comprised  of  zeros.  Rows  two,  four,  amd  six  are  linearly  dependent: 


-2tt 

7 


(Rowj)  +  ^(Row4)  +  Rowe  =  0 


These  rows  are  comprised  of  the  contrcd  terms  in  the  relative  translational  equations  of 
motion.  This  rdationship  among  rows  two,  four,  amd  six  holds  for  any  orbital  inclination. 
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4.1-1  Eigenvalues  and  Stability  The  eigenvalues  of  the  F  matrix  are  the  system 
poles.  They  characterize  the  natural,  unforced  behavior  of  the  system.  A  pole  with  an 
imaginary  component  produces  oscillatory  behavior.  A  pole  with  a  positive  real  part 
causes  system  instability.  We  define  instability  to  mean  that  if  system  is  displaced  from  an 
equilibrium  state  (defined  in  next  section),  the  system  will  not  return  to  that  equilibrium 
state  and  remain.  Mathematica  (12)  is  used  to  determine  the  eigenvalues  (A<)  of  F  (twelve 
total): 


A.  =  0,0, 

±  jn,  ijn 

±  n[3A^  -  2AB  -  2AC  +  BC  - 

±  [OA'*  -  12A*B  -  12A^B^  -  12A^C  +  SOA^BC  +  12AB^C 
-  2A^C*  -  16ABC^  +  B^C^  +  4AC®  -  2BC^  +  /2AB]^ 

We  note  here  that  absence  or  presence  of  three  particular  entries  of  the  F  system 
matrix  has  no  effect  upon  the  eigenvalues  of  F.  These  terms  are  —Zn^z/ro,  — 3ni/2ro  and 
—Sn^yfrQ.  These  are  the  so-called  “extra”terms  whose  origin  and  effect  upon  the  equations 
of  motion  were  discussed  in  section  2.2.2.I. 

The  repeated  poles  at  zero  are  the  source  of  the  constant  and  lineair  terms  in  the 
solution  to  the  coupled  x  and  y  translational  equations  of  motion.  They  medre  the  system 
unstable  in  x  and  y.  The  four  purely  imaginary  poles  lead  to  oscillatory  behavior.  One 
pair  is  associated  with  the  coupled  x  and  y  equations  and  the  other  pair  is  associated 
with  the  z  equation.  The  p^es  ±n(3(A  -  B)fC]^  cause  tps  oscillatory  behavior  if  B  >  A. 
Otherwise  the  V'a  equation  is  exponentially  unstable.  The  last  four  complicated-looking 
poles  are  associated  with  the  coupled  tpi,  equations.  They  cause  oscillatory  behavior 
when  C  >  B  >  A.  Otherwise  these  poles  seem  to  cause  instability.  The  POTV  should 
always  be  designed  to  have  C  >  B  >  A  for  patssive  attitude  stability.  There  may  exist 
combinatimis  of  A,  B,  C  where  the  condition  C  >  B  >  A  does  not  hold  where  the  resulting 
poles  are  stable.  This  possibility  was  not  investigated. 
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4.1.2  Equilibrium  Points  System  equilibrium  points  are  values  of  the  state  vector  x 
which  make  i  —  Fx  =  0.  The  solutions  of  Fx  =  0  lie  in  the  nullspace  of  the  matrix  F. 
One  equilibrium  point  is  the  trivial  solution  =  0.  We  note  that  the  nullspace  of  F  can 
not  solely  consist  of  the  zero  vector  since  rows  one  and  four  of  F  are  linearly  dependent. 
Thus  the  rajik  of  the  12  X  12  matrix  F  is  eleven.  The  Mathematica  “NuUSpace”  command 
is  used  to  calculate  the  resulting  equilibrium  point: 


tT 


00  To  000000100 


where  c  is  any  constant.  We  expect  this  equilibrium  point  to  exist  because  in  chapter  two 
we  saw  that  a  tracking  vehicle  ahead  or  behind  the  target  vehicle  in  the  same  orbit  would 
stay  at  that  relative  position  without  any  expenditure  of  control  effort.  The  pitch  angle 
V’a  is  equal  to  y/fo  at  this  equilibrium  point. 


4.1.2  Controllability  Thus  far  we  have  explored  the  unforced  behavior  of  the  sys¬ 
tem.  The  unforced  system  has  been  shown  to  be  unstable.  Control  effort  must  be  expended 
to  maintain  the  system  in  a  non-equilibrium  state.  This  section  addresses  the  issue  of  con¬ 
trollability.  We  find  that  Design  #  1  is  actually  uncontrollable  at  certain  times  in  an 
equatorial  orbit.  Design  #  2  is  found  to  be  fully  controllable  at  all  times  in  an  equatorial 
orbit. 

A  system  is  said  to  be  completely  controllable  if  in  a  finite  time  span  (t  —  to),  any 
desired  system  state  x(t)  can  be  reached  through  some  control  history  M(t)  (3:384-385). 
The  assessment  of  the  controllability  of  a  linear  time-invarismt  system  is  a  simple  process. 
Unfortunately  our  system  is  linear  and  time-varying.  We  will  make  a  simplification  to  the 
system  that  will  allow  us  to  make  a  time-invariant  controllability  analysis  of  both  designs 
in  an  equatorial  orbit. 

How  is  the  controllability  of  a  linear  time-invariant  system  determined?  Begin  by 
forming  the  controllability  matrix  Me'. 


M,  =  [G1FG|---|F-‘G] 
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If  the  rank  of  is  equal  to  the  system  order  JV,  then  the  system  is  completely  controllable. 
Incidentally,  the  observability  matrix  can  be  formed  in  an  analogous  fashion  but  we  do 
not  concern  ourselves  with  this  because  the  combination  of  the  system  matrix  F  and  the 
chosen  output  matrix  H  always  assure  complete  observability. 

Mathematica  can  be  used  to  evaluate  the  rank  of  Mg  when  F  and  G  are  constant 
matrices.  Our  F  is  a  constant  matrix.  The  non-constant  input  matrices  Gi  and  G2  for  the 
two  designs  contain  the  time- varying  expressions  a,  /3,  amd  7.  The  definitions  of  a,  /?,  7 
are  repeated  here  for  convenience: 

a  =  cos^cosucos(fi  —  —  A)  —  cost cos^sin«sin(f2 - 

(3  =  —  cos^sinucos(n  —  —  A)  —  cost  cos^  cos«sin(fi 

7  =  sinicos^  sin(ft  -  —  A) -|- costsin^ 

where  the  time  dependence  of  these  expressions  comes  from  the  quantities  9g  and  u  (see 
eqs.  3.4  and  2.1): 


—  A)  -|-  sin  t  sin  S  sin  u 
-  —  A)  -h  sin  t  sin  S  cos  u 


9g  —  -I-  w©(<  —  to) 

«  =  Uo  -1-  n(t  -  to) 

We  note  that  the  quantities  a,  /?  are  composed  of  trigonometric  functions  with  frequencies 
n  and  a;©.  The  orbital  frequency  n  is  approximately  fifteen  times  larger  than  the  earth’s 
rotational  frequency  a;©  for  a  400  km  orbit. 

We  use  the  “gain-scheduling”  concept  proposed  by  Ridgely  (17)  to  perform  a  “time- 
invariant”  analysis  of  the  controllability  of  our  system.  This  method  “freezes”  a  time- 
varying  system  at  an  instant  in  time,  then  a  time-invariant  controllability  analysis  is  per¬ 
formed  for  that  time  instant.  If  the  system  is  completdy  controllable  at  all  time  instants, 
then  the  time-varying  system  is  completely  controllable.  The  method  becomes  tractable 
when  the  time- varying  system  is  periodic.  Then  the  values  of  all  entries  in  the  input  matrix 
G  are  known  at  any  future  time  t.  Fortunately  our  system  is  periodic  for  any  inclination 
angle.  This  can  be  seen  from  the  periodic  plots  of  the  magnetic  field  in  the  orbital  frame 
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(appendix  A).  The  gain-scheduling  method  is  used  in  this  section  as  weU  as  in  section  4.2 
where  it  is  used  to  discuss  the  design  of  a  time- varying  controller. 

The  periodicity  of  the  magnetic  field  in  orbital  frame  coordinates  is  conceptually 
simplest  when  the  tracking  vehicle  is  traveling  in  an  equatorial  orbit.  Then  the  in-plane 
Cl  amd  Ci  components  of  the  B_  field  are  single-frequency  sinusoids  and  the  out-of-plane 
C3  component  is  constant  (see  appendix  A).  The  orbital  inclination  angle  i  is  zero  and 
the  right  ascension  of  the  ascending  node  Q  is  undefined  for  an  equatorial  orbit.  For  our 
purposes,  we  pick  Q  to  be  zero.  After  the  application  of  the  trigonometric  identity  (eq. 
3.6),  the  expressions  a,  /?,  7  for  the  equatorial  orbit  become: 

Oq  =  cos^cos(«  —  —  A) 

/3o  =  —  cos  S  sin  («  —  @,  —  A) 

7o  =  sin  S 
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Design  #  2  G2(t): 
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L  2^  2^  iSy  «  U  U  ”  J 

(4.6) 

where  the  abbreviations  c  =  cos ,  «  =  sin ,  and  9  =  (u  —  9f  —  X)  have  been  used  for  brevity. 

The  controllability  of  these  conductor-only  designs  in  non-equatorial  orbits  is  not 
mathematically  investigated  in  this  thesis.  However,  we  expect  that  the  designs  have 
contrcdlability  characteristics  similar  to  those  possessed  in  equatorial  orbits  fur  orbital 
inclinations  other  than  11.5*  which  are  less  than  78.5*.  In  appendix  A  we  find  that  at 
certain  times  (once  per  sidereal  day)  the  in-plane  component  (^|  =  BiCi  -I-  B2C2)  of  H 
approadies  zero  as  the  orbital  inclination  approaches  11.5*.  At  t  =  11.5*,  B^  disappears 
entirely  once  per  sidereal  day  .  Thus  in  an  11.5*  indination  orbit  all  the  control  terms 
involving  a  smd  /3  in  the  equations  of  motion  will  disappear  once  per  sidereal  day.  11.5*  is 
the  approximate  angle  between  the  magnetic  p<des  amd  the  corresponding  geographic  poles. 
Also  in  appendix  A  we  And  that  the  out-of-piane  component  {B_^^)  of  B  is  zero  once  per 
sidereal  day  at  orbital  inclinations  greater  than  78.5*.  Thus  all  contrcd  terms  involving  7 
in  the  equations  of  motion  will  disappear  once  per  sidereal  day  at  inclinations  greater  than 
78.5*.  Later  we  And  that  some  of  the  restrictions  mentioned  in  this  paragraph  disappear 
when  a  pair  of  electromagnetic  thrusters  is  added. 


4-9 


Because  both  of  these  designs  are  periodic  in  nature,  we  do  not  have  to  assess  the 
equatorial-orbit  controllability  of  the  designs  at  every  time  instant  in  order  to  prove  con¬ 
trollability  for  all  time.  Instead  we  evaluate  the  rank  of  the  controllability  matrix  for 
both  designs  for  the  cases: 

•  9  =  {u  —  Bg  —  \)  =  2fc(90‘’),  fc  =  0, 1, 2,  ...=>•  sin  0  terms  in  G{i)  go  to  zero 

•  9  =  (u  —  9g  -  X)  =  (2k  -b  1)(90*),  A;  =  0, 1, 2, . . .  =>  cos  9  terms  in  G(t)  go  to  zero 

•  5  =  («  —  -  A)  =  all  other  values  =>■  cos  sin  0  terms  in  G(t)  are  non- zero 

Evaluating  the  equatorial-orbit  controllability  for  these  three  cases  is  equivalent  to  investi¬ 
gating  the  controllability  for  all  time  instants.  We  want  a  design  which  proves  to  be  totally 
controllable  for  all  values  of  (u  —  9g  —  X).  We  use  the  CONTROLLER.M  program  to  eval¬ 
uate  controllability  by  determining  the  rank  of  the  controllability  matrix  Me-  Design  ^  1 
becomes  uncontrollable  when  —  A)  is  an  odd  multiple  of  90®.  Design  #  2  is  completely 

controllable  for  all  valxies  of  (u  —  9g  —  A).  This  is  an  important  result.  We  foresaw  a  result 
like  this  in  chapter  three  when  we  noted  that  the  V*!,  ^2  equations  for  Design  #  1  (eqs.  3.15 
and  3.16)  did  not  include  any  “constant”  control  terms.  Therefore  all  further  discussion 
will  be  limited  to  Design  #  2.  Although  we  have  only  proved  Design  #  2  controllability 
for  an  equatorial  orbit,  we  expect  this  result  applies  to  other  low-inclination  orbits  as  well. 

Surprisingly,  Design  #  2  turns  out  to  be  completely  controllable  at  all  points  in  an 
equatorial  orbit  even  when  some  conductors  are  removed  (columns  of  G  eliminated).  We 
believe  that  any  three  conductors  can  be  removed  while  maintaining  complete  control¬ 
lability.  This  conclusion  arises  from  much  experimentation  with  the  CONTROLLER.M 
program.  This  redundancy  can  increase  reliability  and  free  conductor  towers  for  other  uses. 
There  are  practical  constraints  on  which  conductors  could  be  removed.  For  instance,  use 
of  out-of-plane  cross  conductors  for  gross  translation  is  not  recommended  due  to  structural 
design  limitations. 

4.2  Controller  Design  and  Performance 

In  this  section  we  use  ..he  gain-scheduling  concept  to  design  a  control  system  that 
applies  to  a  particular  instant  in  time.  A  series  of  such  time-specific  controllers  designed 
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Figure  4.1.  Closed-Loop  System  Block  Diagram 


over  the  time  span  of  one  sideriai  day  could  be  combined  to  form  a  full-blown  time- varying 
control  system. 

The  control  system  should  adlow  the  tracking  vehicle  to  stay  a  fixed  distance  away 
from  and  maintain  a  constant  attitude  with  respect  to  the  tairget  vehicle. 

This  goal  is  stated  in  control  terminology  as:  the  tracking  vehicle  must  track  a 
constant  reference  command.  We  use  the  words  “tracking”  and  “command-following” 
interchangeably.  The  reference  command  consists  of  the  variables  we  want  to  track:  3 
relative  positions  and  three  relative  attitude  angles: 

^det  Vdet  V’S*,, 

The  reference  command  vector  r  defined  here  should  not  be  confused  with  the  POTV 
inertial  position  vector  r  defined  in  chapter  two.  The  controller  should  ensure  system 
stability  and  be  robust  with  respect  to  modeling  error  and  noise.  Figure  4.1  shows  the 
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block  diagrajn  of  the  closed-loop  control  system  proposed  to  accomplish  these  goals.  Note 
that  all  states  are  assumed  to  available  for  feedback.  If  this  were  not  the  case,  a  Kalman 
filter  (11:224)  could  be  used  to  obtain  an  estimate  of  the  state.  The  feedback  controller 
K  is  chosen  to  provide  robust  closed-loop  stability.  The  pre-filter  M  is  designed  for  good 
reference  command  tracking.  The  dimensions  of  all  vectors  and  matrices  in  figure  4.1  are 
as  follows: 


^6x1 

= 

reference  command  vector 

^8X1 

= 

control  vector 

^12X1 

= 

system  state  vector 

^6X1 

= 

output  vector 

jpl2xl2 

= 

system  matrix 

qUxO 

= 

input  matrix 

^6x12 

= 

output  matrix 

= 

pre-filter  matrix 

^8x12 

= 

controller  gain  matrix 

There  are  many  methods  of  feedback  controller  design.  We  choose  the  Linear  Quaxlratic 
Regulator  (LQR)  method  because  this  method  guarantees  closed-loop  stability  and  pro¬ 
vides  good  stability  margins  (11:223,227-229).  The  LQR  controller  is  the  controller  K 
that  minimizes  the  performance  index: 


=  /■ 


{^Qx  -L  i^Ru)dt 


subject  to  the  constraint: 

i  =  Fx  -1-  Gu 

where  Q  is  the  state-wdghting  matrix  and  R  is  the  control-weighting  matrix.  These  ma¬ 
trices  are  of  the  form  pU,  where  U  is  the  identity  matrix,  p,  is  the  weighting  on  the  Q 
matrix.  Pu  is  the  weighting  on  the  R  matrix.  These  weightings  are  chosen  by  the  designer. 
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Usually  one  of  the  weightings  is  fixed  at  unity  while  the  other  is  varied.  Then  the  resulting 
output  response  ajid  corresponding  control  behavior  are  studied.  The  methodology  behind 
the  derivation  of  the  controller  K  that  minimizes  the  performance  index  J  is  described  by 
Maciejowski  (11:222-227)  and  will  not  be  gone  into  here.  The  resulting  K  is  a  constant 
gain  matrix.  In  practice,  the  MATLAB  command  “LQR(F,G,Q,R)”  is  used  to  derive  the 
gain  matrix  K. 

After  the  controller  K  is  designed,  we  derive  the  pre-filter  M.  M  should  be  designed 
with  optimal  command-following  in  mind.  We  assume  that  M,  like  K,  is  a  constant  gain 
matrix.  After  studying  the  interconnections  of  the  system  blocks  in  figure  4.1  we  write: 

X  =  Fx  Gu 

=  Fx  +  G(Mr-  Kx) 

=  (F-GK)x  +  GMr 
y  =  Hx 

X  should  be  zero  in  the  steady-state  as  a  consequence  of  good  tracking  of  the  constant 
reference  command.  Then  the  state  equations  become: 

i  =  0  =  {F-GK)x„  +  GMr 

h.  = 

These  equations  combine  to  yield: 

=  -H(F  -  GKy^GMr  (4.7) 

But  since  the  output  should  track  the  reference  command  r,  the  pre-filter  M  must  be 
chosen  to  satisfy: 

-  H{F  -  GKY'^GM  =  U  (4.8) 
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where  U  is  the  identity  matrix.  We  note  that  M  is  non-square.  Solving  this  equation  for 
M  we  arrive  at: 

M  =  -{H(F  -  GK)-^G)^  (4.9) 

where  the  symbol  denotes  the  pseudo-inverse  operation. 

4-S.l  Stability  Robustness  Uncertainty  in  the  plant  model  can  lead  to  closed-loop 
system  instability.  The  quality  of  a  system  that  makes  it  stable  in  the  face  of  uncertainty 
is  called  stability  robustness.  There  is  a  major  source  of  uncertainty  in  the  POTV  plant 
model.  This  is  the  magnetic  field  magnitude  error  inherent  in  the  use  of  the  simple  dipole 
modd  of  the  geomagnetic  field.  This  error  may  be  up  to  30%  at  the  earth’s  surface 
and  decreases  as  orbital  altitude  increases.  In  control  terminology  this  error  is  termed 
a  “multiplicative  perturbation  at  the  plant  input”.  Fortunately  the  LQR  design  method 
guarantees  a  gain  margin  of  6  dB  when  the  weighting  matrices  used  are  of  the  form  pU 
(11:229).  This  is  more  than  enough  margin  to  guarantee  POTV  closed  loop  stability  in  the 
face  of  magnetic  field  modeling  error.  This  margin  also  helps  overcome  model  linearization 
error. 
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^.2.2  Steady-State  Controllability  The  most  surprising  result  from  experimentation 
with  the  CONTROLLER.M  program  is  that  even  though  the  system  is  completely  con¬ 
trollable,  it  does  not  appeau"  to  be  steady-state  controllable  (a  term  coined  for  this  thesis  by 
Dr.  B.  Liebst).  In  other  words,  the  POTV  can  reach  any  position  state  but  can’t  remain 
at  that  state  unless  it  happens  to  be  an  equilibrium  state.  No  controller  can  be  designed 
that  acceptably  tracks  a  reference  command  consisting  of  three  relative  positions  and  three 
attitude  angles  for  a  POTV  using  only  electrical  conductors  for  thrust.  All  attitude  angles 
are  tracked  but  no  position  components  are  tracked  accurately.  The  designed  controller  for 
the  conductor-only  POTV  works  acceptably  only  when  the  CONTROLLER.M  program 
is  modified  to  include  only  two  of  the  position  components  plus  the  three  attitude  angles 
in  the  reference  command  vector.  The  position  component  that  was  not  included  in  the 
reference  command  then  behaves  in  an  unpredictable  manner.  The  problem  arises  from 
the  rank  deficiency  of  the  G  input  matrix  as  shown  below.  Remember  we  have  shown 
already  that  its  rank  can  never  exceed  five  because  rows  two,  four,  and  six  are  linearly 
dependent.  The  rank  deficiency  arises  from  the  “translational  part”  of  the  input  matrix. 
To  understand  more  about  why  G  should  have  rank  six  for  steady-state  controllability,  we 
return  to  equations  4.7  and  4.8: 


y^^  =  -H{F-GK)-^GMr 
-H{F-GK)-^GM  =  U 

Recall  that  the  second  equation  is  the  necessary  condition  for  the  steady-state  output 
to  track  the  six-component  reference  command  r.  The  6x6  identity  matrix  U  obviously 
has  rank  six.  According  to  Strang  (20:201),  the  rank  of  the  expression  -H{F-GK)~^GM 
must  be  less  th’"  or  equal  to  the  smallest  matrix  rank  found  in  the  expression.  We  have 
already  shown  mat  the  rank  of  G  r.  always  less  than  or  equal  to  five.  Thus  the  expression 
-H{F -  GK)~^GM  can  r.e«erhave  rank  six  and  thus  cam  newer  equal  the  identity  matrix. 
The  steady-state  output  y^^  can  only  track  five  of  the  components  of  the  six-component 
reference  command  r. 
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The  seeming  solution  to  the  steady-state  position  controllability  problem  is  to  some¬ 
how  increase  the  rank  of  G  to  six.  We  introduced  conductors  to  the  design  at  many 
different  locations  and  orientations  in  an  attempt  to  accomplish  this  desirable  result.  The 
new  conductors  added  columns  to  the  G  matrix  but  unfortunately  G  still  lost  rank  in  the 
old,  familiar  fashion: 

— ^(Row2)  -b  ^(Row^)  +  Rowg  =  0 

We  also  deleted  conductors  to  make  the  design  unsymmetrical  with  similar  lack  of  success. 


Steady-state  controllability  cam  also  be  interpreted  as  the  ability  of  the  control  system 
to  exert  thrust  in  any  direction  at  any  time.  Clearly  this  is  not  the  case  with  the  POTV. 
The  cross-product  nature  of  the  electrodynamic  force  tells  us  that  a  force  can  not  be 
exerted  in  the  direction  of  the  instantaneous  B  vector: 


Fb=ILxB  (4.10) 

We  thus  make  the  important  conclusion  that  no  conductor-only  design  configuration  Icawis 
to  steady-state  position  controllability.  Thus  while  the  system  is  completely  controllable 
using  only  electrodynamic  forces  (any  state  can  be  reached  in  a  finite  amount  of  time),  it 
cannot  mmntain  any  state  continuously  other  than  equilibrium  states. 

One  method  of  improving  the  rank  of  the  G  matrix  is  to  introduce  force-producing 
devices  which  do  not  rely  upon  the  electrodynamiic  effect.  We  introduce  a  pair  of  electro¬ 
magnetic  thrusters  to  produce  out-of-plane  thrust.  They  are  attached  to  the  centers  of  the 
out-of-plane  crosses  and  produce  thrust  in  the  “north”  or  “south”  out-of-plane  directions. 
Only  one  thruster  operates  at  a  time.  The  modified  control  vector  u  is: 

^  ~  [  ^1.  h.  h,  T 
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tP 


for  a  total  of  niae  controls.  Positive  T  means  that  the  southern  cross  thruster  is  firing  and 
negative  T  mesuis  the  northern  cross  thruster  is  firing.  The  modified  12x9  G{t)  matrix  is: 
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which  has  a  rank  of  six.  The  addition  of  electromagnetic  thrusters  makes  the  POTV 
completely  steady-state  controllable.  In  addition  to  this  desirable  result,  the  addition  of  the 
thrusters  has  another  benefit:  from  experimentation  with  the  CONTROLLER.M  program 
we  expect  that  the  POTV  becomes  completely  controllable  at  all  points  in  an  ll.S"  orbit 
(the  conductor-only  design  was  already  completdy  controllable  for  other  low-inclination 
orbits).  This  was  not  the  case  for  a  conductor-only  design.  Unfortunately,  we  do  not 
bdieve  that  the  thrusters  make  the  POTV  completely  controllable  at  all  points  in  orbits 
with  inclinations  of  greater  than  78.5^  (as  was  also  the  case  with  the  conductor-only  design). 
The  Design  #1  configuration  may  be  a  more  appropriate  design  for  bugh-incUnation  orbits. 
This  possilniity  is  not  investigated. 
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V.  Steady-State  Tracking  and  Docking  Performance 


In  the  previous  chapter  we  discussed  the  design  of  a  controller  that  enables  the  POTV 
to  track  a  target  vehicle.  In  this  chapter  we  discuss  POTV  conductor  sizing,  steady-state 
tracking  performance,  and  docking  performance.  General  features  of  tracking/docking  per¬ 
formance  are  related  to  controller  design  parameters.  One  steady-state  tracking  example 
and  two  docking  examples  are  shown  and  results  described. 


5.1  Vertical  In-Plane  Conductor  Sizing  for  Orbit  Transfer 

In  this  section  we  size  the  vertical  (*)  ia-plane  conductors  for  an  orbital  transfer 
requirement.  The  horizontal  (y)  in-plane  conductors  are  discussed  in  the  docking  section 
of  this  chapter.  The  out-of-plane  cross  conductors  and  their  non-conducting  support  towers 
ayre  not  explicitly  sized.  We  assume  that  the  cross  conductors  are  much  smaller  than  the 
in-plane  conductors  because  we  expect  they  will  be  used  primarily  for  attitude  control. 

A  practical  method  to  move  a  low  thrust  vehicle  like  the  POTV  from  a  low  prograde 
orbit  to  a  higher  altitude  prograde  orbit  is  to  exert  a  small  constant  thrust  along  the 
orbital  velocity  vector  (eastward).  Thrust  in  the  westward  direction  would  cause  the  orbit 
to  degrade.  The  trajectory  resulting  from  a  small  constaint  eastward  thrust  is  a  gentle  spiral 
outward  from  the  lower  orbit  to  the  higher  orbit  assuming  the  lower  orbit  was  circular. 
Wiesel  (24:89-90)  derives  the  time-of-flight  equation  for  this  method  of  orbital  transfer: 


,  _  Mty/GM^  f  1  1 


(5.1) 


where  F  is  the  applied  thrust,  Mt  is  the  total  system  mass,  is  the  lower  orbital  radius, 
and  r„  is  the  upper  orbital  radius.  We  use  this  equation  to  determine  the  thrust  required 
to  move  the  POTV  from  one  circular  orbit  to  another  for  a  given  transfer  time  and  vehicle 
mass. 

The  ASSET  version  of  the  POTV  spacecraft  is  required  to  move  aluminum  salvaged 
from  space  shuttle  external  fuel  tamks  from  an  orbital  altitude  of  approximately  300  km 
to  an  altitude  of  460  km.  Here  we  assume  that  the  space  shuttle  can  deliver  its  external 
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tank  to  300  km  (4).  The  proposed  maximum  altitude  of  the  space  station  Freedom  is 
460  km  (13).  The  space  station  travels  in  an  almost  circular  orbit  (e  =  0.009).  We 
assume  for  ezise  of  derivation  that  the  POTV  is  traveling  in  circular  orbit  when  it  finishes 
scavenging  a  shuttle  external  tank  for  aluminum.  Now  we  can  determine  the  amount  of 
thrust  required  to  move  a  50,000  kg  POTV  from  300  km  to  460  km  in  thirty  days.  We 
have  chosen  50,000  kg  as  the  maximum  allowed  value  of  the  total  system  mass.  This  figure 
includes  a  fixed  31,300  kg  central  mass  and  the  total  mass  of  the  conductor  towers,  which 
has  not  been  determined  yet.  Evaluating  equation  5.1  we  find  the  required  force  (thrust) 
is  approximately  1.75  newtons.  Now  we  will  “size”  the  two  vertical  in-plane  conductors 
(one  “up”  and  one  “down”)  required  to  produce  this  thrust  assuming  a  certain  current 
level  can  be  sustained.  Penzo  and  Ammann  (15:129)  discuss  a  proposed  electrodynamic 
tether  system  capable  of  handling  125  amps  of  current.  We  take  this  to  be  a  reasonable 
sustainable  current  limit.  Then  the  familiar  Lorentz  equation  (eq.  1.1)  is  used  to  make  an 
engineering  estimate  of  the  total  conductor  length  required: 

Fb=ILxB 

We  assume  that  the  needed  force  (1.75  newtons)  is  produced  by  the  interaction  of  the 
two  vertical  in-plane  conductors  with  the  out-of-plane  component  of  the  geomagnetic  field 
{B±).  At  the  space  station  inclination  of  28.5®  (13),  B±  consists  of  a  large  constant 
component  and  a  smaller  sinusoidal  component  (see  appendix  A  for  a  plot  of  Bj_  (B3)  at 
i  =  28.5*).  Over  the  altitude  range  of  300-460  km,  B±  varies  from  a  low  of  approximately 
1.94  •  10“*  tesla  to  a  high  of  2.60  •  10“®  tesla.  For  calculation  purposes  we  take  B±  to  be 
approximately  a  constant  2.26  ■  10~’’  tesla  over  the  orbital  transfer.  Solving  the  Lorentz 
equation  for  L,  we  find  that  the  two  vertical  in-plane  conductors  must  have  a  combined 
length  of  approximately  6S0  meters  for  the  mission  parameters  we  have  outlined. 

We  would  like  to  minimize  the  amount  of  thrust  needed  for  planar  orbital  radial 
change  because  this  minimizes  the  product  of  current  and  conductor  length.  Since  power 
is  proportional  to  current  squared,  current  is  a  variable  worth  minimizing.  There  are 
several  ways  to  accomplish  this: 
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•  Minimize  the  amount  of  radial  distance  to  be  traversed. 

•  Increase  the  amount  of  time  allowed  for  the  maneuver. 

•  Decrease  the  system  mass  by  minimizing  the  central  mass  and  by  minimizing  the 
linear  mass  density  of  the  conductor  towers. 

We  can  also  minimize  the  amount  of  current  needed  for  a  given  thrust  level  by  making  the 
conductors  as  long  as  possible.  Unfortunately  this  conflicts  with  the  goal  of  minimizing 
total  system  mass. 

What  size  power  plant  would  be  needed  for  the  particular  orbital  transfer  example 
(ASSET  POTV)  we  have  described?  We  use  the  approximate  induced  voltage  equation 
(eq.  C.l)  for  low-inclination  orbits,  and  the  power  equation  (eq.  C.2)  derived  in  appendix 
C  to  derive  to  the  necessary  power  supply  power: 

~  101  volts 

where  we  have  used  r  =  300  km  to  be  conservative  in  our  calculation  of  V<„^  because  the 
induced  voltage  is  greatest  at  the  lower  orbital  radius.  Then  we  derive  the  needed  power: 


Pp.=IVi„,  +  PR,,, 


where: 


R 


tot 


P-C  +  Rim  +  R 


load 


ff  we  assume  that  all  power  is  being  used  for  thrust,  Ru,^  is  zero.  Erom  chapter  one  we 
know  Rign  varies  from  approximately  one  ohm  to  twenty  ohms.  The  conductor  resistance 
Rf  is  equal  to  the  conductor  resistance  per  unit  length  multiplied  by  the  total  conductor 
length.  The  125  unp  exsunple  electrodynamic  tether  discussed  by  Penzo  and  Ammann 
has  a  resistance  per  length  of  approximatdy  7.7  - 10~*  ohms  per  meter  (15:129).  Using 
this  value,  the  conductor  resistance  is  0.05  ohms.  So  most  of  the  resistance  arises  in  the 
ionospheric  current  path.  The  total  power  for  this  example  for  a  current  of  125  amps 
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and  the  worst  case  ionospheric  resistance  of  twenty  ohms  is  calculated  to  be  about  325,000 
watts,  with  most  of  the  power  being  dissipated  in  the  ionosphere  resistance.  If  we  use  the 
best  case  one  ohm  value  for  the  ionospheric  resistance,  the  total  power  is  approximately 
44,000  watts.  These  figures  have  a  great  deal  of  built-in  uncertainty  because  the  resistance 
of  the  ionosphere  is  not  well  understood.  To  minimize  power  requirements  a  practical 
vehicle  must  be  designed  to  minimize  all  sources  of  resistance.  More  importantly,  conductor 
lengths  should  be  maximized  to  reduce  the  amount  of  current  needed  for  a  given  thrust 
level  since  power  is  proportional  to  the  current  squared.  A  possible  option  to  maximize 
conductor  length  would  be  to  deploy  retractable,  low-mass,  flexible  tethers  at  the  end  of 
the  rigid  conductors.  These  flexible  tethers  would  be  suitable  for  the  orbital  transfer  task 
since  precision  position  control  is  not  required. 

Would  the  magnetic  field  generated  by  a  125  amp  current  in  the  vertical  in-plane 
conductors  interfere  with  the  currents  flowing  through  other  conductors?  If  we  assume 
a  minimum  distance  between  conductors,  we  can  answer  this  question.  Let’s  assume  a 
minimum  separation  distance  of  5  meters.  Then  equation  1.2  is  used  to  estimate  the 
magnetic  field  generated  by  a  125  amp  current  at  a  5  meter  radial  distance.  The  answer  is 
5- 10“®  tesla.  This  is  approximately  one-fourth  of  the  geomagnetic  strength  at  the  altitudes 
at  which  we  expect  to  use  the  ASSET  POTV.  A  robust  control  system  should  be  able  to 
handle  this  “disturbance”. 

5.2  Steady-State  Tracking  Performance 

The  controUer  design  progrsun  CONTROLLER.M  is  a  MATLAB  program  written 
for  this  study  which  lets  the  user  enter  various  parameters  which  describe  the  POTV 
system,  enter  the  initial  standofi  position  and  relative  attitude  (initial  state),  enter  the 
desired  standoff  position  and  rdative  attitude  (reference  commamd),  enter  the  LQR  design 
weightings  and  Pu,  and  evaluate  the  resulting  transient  and  steady-state  response  of 
the  outputs  and  controls.  See  appendices  E  and  F  for  three  sample  runs,  corresponding 
plots,  and  a  complete  program  listing.  All  conclusions  reached  about  steady-state  and 
non-steady-state  tracking  performauice  in  this  and  following  sections  are  the  result  of  ex¬ 
perimentation  with  this  program.  The  input  matrix  G{t)  in  the  program  includes  the  pair 
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of  electromagnetic  plasma  thrusters  added  at  the  end  of  chapter  four  as  well  as  the  Design 
^  2  conductor  configuration.  We  evaluate  the  tracking  performance  of  the  controller  for 
two  general  cases: 

•  Steady-State  Tracking:  r  =  x(to)  The  reference  command  is  equal  to  the  initial  state 
of  the  system.  The  POTV  is  being  directed  to  stay  at  its  initial  position  and  attitude 
relative  to  the  target  vehicle. 

•  Non-Steady-State  Tracking:  r  ^  *(to)  The  reference  command  is  not  equal  to  the 
initial  state  of  the  system.  The  system  is  being  commanded  to  a  new  relative  position 
and/or  rdative  attitude.  This  is  the  case  when  the  POTV  is  commanded  to  dock 
with  a  taurget  vehicle.  We  show  two  docking  examples. 

In  this  section  we  discuss  the  steady-state  tracking  case.  The  docking  cases  are  discussed 
in  the  next  section. 

A  steady-state  tracking  example  is  the  first  of  three  examples  shown  in  appendix 
E.  A  sample  run  at  the  beginning  of  the  appendix  shows  the  data  the  user  enters  to 
describe  the  POTV  configuration,  design  weightings,  and  tracking  requirements.  Plots 
are  made  of:  position  states,  attitude  states,  forces,  moments,  and  conductor  currents. 
For  the  first  example  we  choose  an  length  of  310  meters  to  provide  the  required  large 
steady-state  force  in  the  x  direction.  This  is  also  the  length  derived  in  the  previous  section 
for  the  orbital  transfer  requirement.  We  have  assumed  here  without  proof  that  the  POTV 
csm  rotad;e  90  de^ees  about  its  z  axis  when  it  transitions  from  orbital  transfer  mode  to 
steady-state  tracking  mode.  From  section  2.2.1. 2  we  know  that  no  steady-state  force  will 
be  required  in  the  y  direction.  Thus,  li  can  be  chosen  to  be  “short”.  We  aissume  here 
that  most  of  the  generated  forces  are  the  result  of  current  interactions  with  the  large  out- 
of-plane  component  of  B_  present  in  low-inclination  orbits.  We  choose  an  li  length  of  50 
meters  because  this  is  a  convenient  length  if  the  vertical  in-plane  towers  are  also  used  for 
docking  purposes.  The  non-conducting  out-of-plane  towers  are  chosen  to  be  50  meters  long 
amd  all  the  conductor  lengths  in  the  out-of-plane  crosses  are  chosen  to  be  20  meters  long. 
The  out-of-plane  crosses  are  designed  with  attitude  contrcfi  in  mind.  We  do  not  want  large 
currents  to  be  generated  in  the  crosses  because  the  resulting  large  forces  might  cause  the 
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non-conducting  support  towers  to  fail  structurally.  All  of  the  conductor  length  information 
presented  here  and  other  information  about  the  specific  vehicle  used  in  the  examples  is 
tabulated  in  appendix  D.  The  same  vehicle  configuration  is  used  for  all  three  examples. 

The  LQR  controller  design  luethodolc^  we  have  used  appears  to  be  well-suited  to 
the  steady-state  tracking  example.  The  plots  in  appendix  section  E.l  backup  the  following 
observations;  No  transient  position  output  behavior  is  observed.  The  position  outputs 
aure  always  constant  and  equal  to  their  initial  values  x(to)  regardless  of  state  and  control 
weighting.  The  attitude  outputs  each  exhibit  a  decaying  oscillatory  sort  of  behavior  which 
settles  down  to  a  constant  steady-state  value.  The  observed  sign  and  magnitude  of  the 
steady-state  forces  and  moments  are  the  same  as  those  predicted  in  sections  2.2. 1.2  and 
2.2.2.2.  If  we  increase  the  state-weighting  the  amplitude  of  oscillation  of  the  attitude 
outputs  decreases  and  settling  times  decrease.  We  believe  that  the  extremely  small  non¬ 
constant  force  observed  in  the  y  direction  arises  from  moment  required  for  attitude  control. 
We  think  that  other  small  non-constant  forces  required  for  attitude  control  do  not  appear 
on  the  plots  of  forces  in  the  x  and  z  directions  because  the  constant  forces  needed  for  steady- 
state  position  standoff  are  many  orders  of  magnitude  greater  than  the  non-constant  forces 
required  for  attitude  control. 

5.3  General  Non-Steady-State  Tracking  Performance 

The  findings  in  this  section  apply  to  both  of  the  docking  examples  that  follow.  Unlike 
steady-state  tracking,  transient  output  behavior  is  observed  for  position  as  well  as  attitude 
outputs.  All  outputs  exhibit  the  ‘‘decaying  oscillation  to  a  constant  steady-state  value” 
type  behavior.  As  in  steady-state  tracking,  the  amplitude  of  oscillation  and  the  settling 
time  decrease  as  the  state- weighting  />.  is  increased.  Unfortunately  as  is  increased, 
the  transient  initial  control  requirement  increases.  We  adso  observe  that  as  the  difference 
between  the  initial  state  vector  and  the  reference  command  vector  increases,  the  transient 
initial  control  requirement  increases.  Also  for  a  fixed  state-weighting,  as  the  difference 
increases,  the  amplitude  of  the  initial  oscillation  of  the  outputs  increases. 

Heavy  state-weighting  is  wanted  for  desirable  transient  response  behavior.  Light 
state-weighting  is  desired  to  minimize  the  transient  initial  control  effort  required.  This  is 
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a  good  example  of  the  conflicting  requirements  commonly  found  in  engineering  problems. 
One  solution  is  to  make  the  state- weighting  p*  large  enough  to  satisfy  the  speed  of  response 
requirements  of  a  particular  problem.  Then  the  initial  reference  command  is  made  equal 
to  the  initial  state.  Then  the  reference  command  is  moved  away  from  the  initial  state 
toward  the  desired  final  state  by  small  enough  increments  such  that  the  initial  control 
requirement  is  brought  below  an  acceptable  level.  We  refer  to  this  process  as  “stepping” 
the  reference  command.  Another  possible  solution  to  this  problem  which  is  applicable  to 
R-BAR  docking  is  presented  in  the  R-BAR  docking  subsection. 

5.4  Docking  Examples 

In  the  docking  examples  in  this  section  we  assume  that  the  POTV  has  been  moved 
close  enough  to  the  target  vehicle  for  the  linearized  relative  positional  and  rotational  equa¬ 
tions  of  motion  to  be  valid.  Examples  are  shown  in  appendix  E  for  two  different  in-plane 
docking  approaches.  The  vehicle  lengths  used  for  the  steady-state  tracking  example  are 
used  here.  An  out-of-plane  docking  approach  (Z-BAR)  is  discussed. 

5.4.1  R-BAR  An  R-BAR  approach  assumes  that  both  of  the  vehicles  are  in  the 
same  orbital  plane.'  The  POTV  approaches  the  target  vehicle  by  moving  along  the  radius 
vector  connecting  the  center  of  the  earth  and  the  target  vehicle.  Forces  in  the  x  and  y 
direction  are  exerted  to  move  the  POTV  downward/upward  along  the  rsulius  vector  and  to 
keep  the  POTV  from  moving  east  or  west  as  it  moves  along  the  radius  vector.  A  specific 
example  is  shown  in  appendix  section  E.2. 

R-BAR  position  stmidoff  in  the  x  direction  is  the  most  energy-expensive  type  of 
standoff.  Standoff  in  the  z  direction  requires  three  times  less  control  effort.  Standoff  in  the 
y  direction  requires  no  control  effort.  We  have  repeated  these  conclusions  here  from  section 
2.2. 1.2  for  convenience.  Long  horizontal  in-plane  conductors  are  required  to  produce  the 
necessary  force  in  the  x  direction.  We  can  in  fact  make  an  engineering  approximation 
to  the  length  of  the  horizontal  in-plane  conductors  (/2)  needed  for  a  given  vehicle  mass, 
maximum  current,  and  R-BAR  standoff  distance  (x-direction).  The  following  derivation 
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makes  use  of  equations  1.1  and  2.17. 


l2{2l2)Bx  =  3n^{Mcyl  +  2/»(/i  +  /j  +  ^*  +  ^e,  + 


which  can  be  solved  for  I2  to  yield; 

7  _  +  ^P(h  +  hi  +  ht)  /c 

^  “  2l2Bj_/Zn^x.,  -2p  ^ 

We  note  that  I2  is  extremely  sensitive  to  the  denominator  expression  2l2Bxl3n?x,,  —  2p. 
Any  parameter  in  this  expression  that  varies  will  significantly  affect  the  sizing  of  /j.  We 
repeat  some  of  the  example  vehicle  descriptive  information  from  appendix  D  here  for  the 
calculation  of  l2-  Using  Mcyi  =  31300  kg,  p  =  7.3  kg/m,  li  =  50  m,  /,  =  50  m,  !<.,  =  Ic,  =  20 
m,  I2  =  125  amperes,  —  2.6  •  10"®  tesla,  n*  =  1.28  •  10"®  rad^/sec*,  and  x,,  =  100  m: 

I2  ^  14000  m 

which  is  clearly  an  excessive  length  for  “rigid”  conductor  towers.  If  we  recalculate  I2  using 
an  X  standoff  distance  of  50  m,  the  required  I2  length  is  approximately  1700  m.  One  way 
the  problem  of  long  I2  conductors  can  be  remedied  is  by  allowing  more  current.  This  could 
be  done  by  “bundling”  several  conductors  together  in  a  single  structural  tower,  splitting 
the  required  current  among  several  conductors. 

5.4-S  V-BAR  When  the  POTV  and  the  target  vehicle  are  in  the  same  orbit,  one 
way  the  POTV  can  approach  the  target  is  by  moving  in  the  y  direction  along  the  orbit. 
This  is  known  as  a  V-BAR  docking  approach.  Forces  in  the  x  and  the  y  direction  are 
exerted  to  move  the  POTV  eastward/westward  in  the  orbit  and  to  keep  the  POTV  from 
moving  up  or  down  in  the  orbital  plane.  A  specific  example  of  V-BAR  docking  is  shown 
in  appendix  section  £.3. 
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5.4-3  Z-BAR  The  Z-BAR  approach  starts  with  the  POTV  and  target  vehicle  not 
in  the  same  orbital  plane.  The  POTV  approaches  the  target  vehicle  along  the  normal  to 
the  target  orbit.  An  example  for  this  approach  is  not  shown  in  this  thesis  because  when 
we  ran  a  Z-BAR  approach  with  CONTROLLER.M,  no  currents  were  used  for  position 
control.  The  only  “controls”  that  were  used  for  position  change  were  the  electromagnetic 
plasma  thrusters.  The  conductor  currents  were  very  small  and  appeared  to  be  used  only 
for  attitude  control. 
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VI.  Conclusions  and  Recommendations 


In  this  final  chapter  we  summarize  the  major  results  of  the  thesis  work  and  make 
recommendations  for  further  research. 

While  the  derivation  of  the  linearized  relative  position  equations  of  motion  (Clohessy- 
Wiltshire  equations)  has  been  “done”,  we  are  not  aware  of  others  performing  our  particular 
derivation  of  the  small-angle  relative  rotational  equations  of  motion.  These  equations  are 
the  rotational  equivalent  of  the  Clohessy- Wiltshire  equations.  From  these  two  sets  of  equa¬ 
tions  we  derived  the  steady-state  forces  and  moments  required  to  maintain  a  continuous 
standoff  in  position  and  attitude.  These  steady-state  predictions  were  confirmed  in  the 
example  force  and  moment  plots  contained  in  appendix  E.  Further  research  should  be 
done  to  determine  at  what  point  the  linearization  assumptions  (used  extensively  in  the 
derivation  of  the  equations  of  motion)  start  to  break  down. 

The  derivation  of  the  geomagnetic  field  in  spacecraft  coordinates  was  the  most  im¬ 
portant  early  result  in  the  thesis  research.  Without  this  knowledge  the  derivation  of  the 
full  equations  of  motion  including  magnetic  forces  and  moments  would  have  been  impossi¬ 
ble.  The  behavior  of  the  in-plane  and  out-of-plane  B_  components  over  the  orbital  period 
and  over  the  period  of  a  day  was  unexpected.  IVrther  research  should  employ  a  better 
model  of  the  geomagnetic  field  than  the  simple  dipole  model  to  determine  if  the  in-plane 
components  of  B_  really  do  “disappear”  at  one-day  intervals  in  an  11.5"  inclination  orbit.  A 
similar  analysis  could  be  done  for  the  out-of-plane  component  of  in  the  78.5"  inclination 
orbit.  The  error  arising  from  our  use  of  the  simple  dipole  model  was  compensated  for  by 
the  design  of  a  robust  control  system.  NASA  software  (15:168)  exists  that  more  accurately 
models  the  geomagnetic  field  in  geocentric  coordinates.  This  software  could  be  modified 
to  translate  between  geocentric  coordinates  and  spacecraft  coordinates  if  this  has  not  been 
done  already. 

The  derivation  of  B,  in  spacecraft  coordinates  enabled  us  to  derive  the  full  equations 
of  motion  including  the  effects  of  magnetic  force  as  well  as  gravitational  force.  While 
aerodynamic  and  oblateness  effects  were  described  in  appendix  B,  they  were  not  considered 
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in  the  controller  design.  These  forces  axe  significant  in  low  earth  orbit  and  are  deserving 
of  further  study.  The  air  drag  especially  needs  to  be  quantified. 

After  deriving  the  full  equations  of  motion  we  used  a  state  space  model  to  represent 
the  system.  This  sort  of  representation  was  convenient  in  addressing  the  issues  of  control¬ 
lability,  steady-state  controllability,  and  controller  design.  We  found  that  all  conductor 
designs  were  not  equal.  Design  #  1  lost  control  authority  in  equatorial  orbits  at  points 
where  Design  ^  2  did  not.  Future  researchers  could  investigate  whether  Design  ^1  is  a 
better  design  than  Design  ^2  for  high-inclination  orbits.  We  suspect  that  this  is  the  case. 
We  also  found  that  some  of  the  conductors  were  redundant.  Further  work  should  address 
the  design  of  a  minimal  conductor  configuration  with  one  or  more  redundant  conductors 
for  docking  purposes. 

One  of  the  most  significant  findings  of  the  thesis  was  the  seeming  impossibility  of 
building  an  electrical  conductor-only  design  that  was  steady-state  controllable.  This  was 
clearly  a  disappointing  result.  We  had  hoped  to  build  a  conductor-only  design.  The 
axldition  of  a  pair  of  electromagnetic  plasma  thrusters  allowed  the  POTV  to  track  all 
components  of  position  as  well  as  attitude. 

An  important  accomplishment  of  this  work  was  the  successful  design  of  a  closed-loop 
controller  using  the  LQR  method.  This  was  not  a  “given”  at  the  beginning  of  the  research 
because  the  open-loop  state  responses  of  the  system  to  the  control  currents  appeared  to  be 
hopelessly  coupled.  In  other  words,  a  control  action  giving  a  desired  change  in  one  state 
would  produce  simultaneous  undesired  changes  in  other  states. 

The  LQR  design  method  had  one  major  shortcoming:  no  bounds  could  be  placed 
on  the  amount  of  control  energy  the  controller  could  “ask  for.”  Future  researchers  could 
conceivably  use  a  design  methodology  other  that  LQR  that  would  allow  the  designer  to 
place  bounds  upon  the  allowable  contrcd  energy.  Also,  the  possibility  of  the  design  of  a 
reduced  state  controller  should  be  investigated. 

We  felt  that  the  gain-scheduling  concept  was  jui  adequate  way  to  use  time-invariant 
control  techniques  to  analyze  a  time-varying  system.  Future  researchers  could  modify  the 
CONTROLLER.M  program  to  produce  a  controller  A'/pre- filter  M  combination  for  many 
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“steps”  (times)  in  the  POTV  orbit  over  the  period  of  a  siderial  day.  Then  each  component 
of  the  gain  matrix  K  could  be  plotted  versus  time  to  deduce  the  characteristics  of  the 
time-varying  controller.  We  predict  that  the  components  of  the  time-varying  gain  matrix 
K{t)  will  consist  of  constants  and  sinusoidal  terms. 

An  important  result  of  this  research  was  the  finding  that  large  tracking  vehicles 
require  large  currents  for  steady-state  position  standoff  in  amy  direction  other  than  along 
the  orbital  velocity  vector.  Perhaps  this  problem  could  be  handled  by  “bundling”  several 
conductors  together  in  a  single  structural  tower,  splitting  the  required  current  among 
several  conductors.  This  idea  may  not  be  technically  feasible.  The  currents  generated  by 
the  PMGs  associated  with  each  conductor  might  interfere  with  each  other.  This  idea  is  a 
good  candidate  for  further  research. 

We  addressed  the  problem  of  orbital  radial  change  in  the  sizing  of  the  vertical  in¬ 
plane  conductors.  The  most  important  type  of  orbital  maneuver  for  the  POTV  after  radial 
change  is  that  of  inclination  change.  Penzo  and  Ammann  claim  that  an  electrodynamic 
tether  can  be  used  to  alter  all  six  orbital  elements  (15:168-169).  Their  ideas  should  be 
applicable  to  the  POTV.  The  amount  of  current  needed  for  maneuvers  decreases  with 
increasing  magnetic  field  strength.  A  Jupiter  Inner  Magnetosphere  Maneuvering  Vehicle 
(15:66-67)  has  been  proposed  to  explore  the  Jovian  system,  where  the  field  strength  is 
approximately  twenty  times  that  of  the  earth.  Another  more  neu  term  application  is 
that  of  pure  attitude  control.  V'e  do  not  know  of  any  research  being  performed  in  the 
area  of  attitude  control  using  multiple,  external  conductors.  A  POTV-like  vehicle  could 
theoretically  be  used  anywhere  a  strong  magnetic  fidd  and  plasma  exist.  We  believe  an 
exciting  application  would  be  to  use  POTV  te<^nology  to  design  a  vehicle  to  explore  the 
magnetosphere  of  the  sun.  On  that  futuristic  note  we  md  the  thesis. 
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Appendix  A.  Plots  of  the  Geomagnetic  Field  in  the  Orbital  Reference  Frame 


These  plots  of  magnetic  field  strength  in  c  frame  coordinates  assume  a  circular  orbit 
of  400  km  altitude.  The  orbital  period  at  this  altitude  is  approximately  5500  seconds. 
The  time  duration  of  the  plots  is  two  siderial  days  (172328  mean  solar  seconds).  Plots  are 
generated  for  orbital  inclinations  of  0®,  5",  11.5®,  28.5®,  60®,  78.5®,  and  90®.  The  inclinations 
of  0®,  11.5®,  28.5®,  and  78.5®  are  especially  important  to  the  thesis  development.  Bi{t) 
and  J?2(0  the  in-plane  components  of  the  magnetic  field.  Bz{t)  is  the  out-of-plane 
component. 
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Appendix  B.  Aerodynamic  and  Earth  Oblateness  Effects  in  Low  Earth  Orbit 


In  chapters  two  and  three  we  quantified  the  gravitational  and  magnetic  forces  and 
torques  acting  on  the  tracking  vehicle.  These  are  the  only  forces/moments  included  in  the 
POTV  equations  of  motion.  Other  forces  and  torques  include  those  caused  by  aerodynamic 
drag,  the  oblate  shape  of  the  earth,  solar  radiation  pressure,  solar  and  lunar  gravity,  dust 
impingement,  etc.  We  believe  the  most  significant  perturnations  in  low  earth  orbit  (LEO) 
are  aerodynamic  drag  and  earth  oblateness. 

The  drag  due  to  friction  with  the  earth’s  atmosphere  opposes  the  orbital  motion  of 
the  POTV.  The  net  effect  of  this  force  is  to  circularize  an  elliptical  orbit  and  to  decrease  the 
orbital  radius.  Wiesel  (24:83),  and  Penzo  and  Ammann  (15:144),  model  the  aerodynamic 
drag  ^  as: 

£Lt  =  -\CiApvl^tC2 

where  Cd  is  the  drag  coefficient,  A  is  the  presented  area  of  the  POTV,  p  is  the  atmospheric 
density,  and  Vrti  is  the  velocity  of  the  POTV  relative  to  the  air.  At  the  orbital  altitude 
range  of  the  POTV  (300-460  km),  a  simple  exponential  model  of  the  atmospheric  density 
p  does  not  suffice.  Penzo  and  Ammann  give  an  empirical  formula  for  the  atmospheric 
density  for  altitudes  greater  than  200  km  when  the  particle  density  is  greater  that  10'^ 
particles  per  cubic  meter: 

1.47  •  10-^«re,(3000  -  T„) 

(1  +  2.9(^4i^))io 

where  p  is  measured  in  kg/meter®,  T  is  the  atmospheric  temperature  in  degrees  Kelvin, 
and  Alt  is  the  orbital  altitude  in  kilometers. 

Aerodynamic  drag  is  the  major  force  to  be  countered  for  simple  station-keeping  in 
low  earth  orbit  (LEO).  No  attempt  is  made  here  to  determine  the  amount  of  vehicle  thrust 
needed  to  counter  the  drag  force.  The  calculation  of  the  drag  coefficient  for  the  complex 
shape  of  the  POTV  is  difficult.  In  practical  tracking  and  docking  applications  the  drag 
coefficient  for  both  the  tracker  and  target  should  be  calculated. 
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The  effects  of  the  oblate  shape  of  the  earth  must  also  be  considered  for  a  vehicle 
in  LEO.  The  earth’s  equatorial  bulge  affects  an  elliptical  orbit  in  two  significant  ways 
(24:86-87): 

•  Regression  of  the  line  of  nodes. 

•  Rotation  of  the  line  of  apsides  (also  known  as  the  advance  of  the  perigee). 


The  second  of  these  effects  does  not  apply  to  a  circular  orbit.  We  have  always  assumed 
that  the  POTV  is  traveling  in  circular  orbit  or  tracking  another  vehicle  in  a  circular  orbit. 
The  regression  of  the  line  of  nodes  is  quantified: 


ZnJ2r% 

2a2(l  -  e2)2 


cost 


which  for  a  circular  orbit  simplifies  to: 


6  — 

2r2 


cost 


where  tl  is  the  time  rate  of  change  of  the  right  ascension  of  the  ascending  node,  n  is  the 
orbital  angular  velocity  (or  mean  motion),  J2  —  0.001082  is  a  number  describing  the  oblate 
shape  of  the  earth,  r®  is  the  radius  of  the  earth,  r  is  the  POTV  orbital  radius,  and  t  is  the 
orbited  inclination.  The  way  this  effect  could  be  included  in  the  POTV  equations  of  motion 
would  be  to  substitute  the  expression  “fio  +  wherever  “fl”  appears  in  the  equations. 


Appendix  C.  Electrodynamic  Propulsion  Power  Requirements 


What  are  the  power  requirements  inherent  in  electrodynamic  propulsion?  The  volt¬ 
age  induced  across  a  conductor  moving  in  a  magnetic  field  is  an  important  variable  in 
power  calculations.  The  voltage  induced  across  a  straight  conductor  moving  through  a 
uniform  magnetic  field  is  given  by  the  relationship  (15:119): 

^ind  =  ^nd  '  Id. 


where  is  the  electrical  field  induced  across  the  conductor  and  v  is  the  velocity  of  the 
conductor  relative  to  the  magnetic  field.  The  velocity  v  of  an  earth  orbiting  conductor 
relative  to  the  geomagnetic  field  is  equal  to  its  orbital  velocity  minus  the  velocity  of  the 
co-rotating  magnetic  fieid.  For  a  circular,  equatorial  orbit,  the  relative  velocity  is  written: 


V  = 


-  w®r 


where  G  is  the  gravitational  constant,  is  the  mass  of  the  earth,  r  is  the  orbital  radius, 
and  Uq  is  the  earth’s  rotational  angular  velocity.  This  expression  is  approximately  correct 
for  low-inclination  orbits.  The  first  term  in  this  equation  is  always  greater  than  the  second 
term  for  low  earth  orbit.  We  note  that  no  voltage  is  induced  if  the  conductor  is  oriented 
such  that  the  induced  electrical  field  =  u  x  ^  is  perpendicular  to  the  length  vector. 
When  the  induced  electrical  field  is  parallel  to  L,  the  induced  voltage  is  simply: 


w^r)BL 


(C.l) 


This  approximate  equation  is  used  in  chapter  five  to  estimate  the  power  requirements  for 
a  radial  orbital  change  application. 

Figure  C.l  shows  an  electrical  circuit  model  of  a  conductor  (tether)  operating  as 
a  thruster.  The  induced  voltage  is  modeled  as  a  voltage  source  with  polarity  opposite 
to  that  of  the  power  supply  voltage  source.  We  note  that  the  current  would  flow  in  the 
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opposite  direction  of  that  shown  if  the  power  supply  voltage  was  less  than  the  induced 
voltage.  This  is  in  fact  the  case  when  a  tether  is  used  for  power  generation.  Unfortunately, 
power  generation  is  always  accompanied  by  electrodynamic  drag  for  a  conductor  in  low 
earth  orbit.  We  are  interested  in  the  power  required  to  operate  a  conductor  as  a  thruster. 
Regarding  figure  C.l  we  use  KirchofF’s  voltage  law  to  write: 

—  ^ind  +  f  Rtot 

where  is  the  power  supply  voltage,  Vin4  is  the  induced  voltage  across  the  conductor, 
/  is  the  current  flowing  through  the  conductor,  and  Rtot  is  the  total  electrical  resistance. 
Rtot  is  the  sum  of  the  resistances: 

Rtot  =  Re  +  R  ion  +  fiload 
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where  R^,  is  the  conductor  resistance,  Ri^n  is  the  ionospheric  resistance,  and  R\oad  is  the 
load  resistance.  Penzo  and  Ammann  state  that  the  ionospheric  resistance  may  be  on  the 
order  of  one  to  twenty  ohms  (15:125).  Actual  orbital  tests  are  needed  to  verify  this  figure. 
The  load  resistance  Rioad  is  derived  from  POTV  equipment  other  than  the  conductors  that 
need  electrical  power.  The  conductor  resistance  Rc  is  equal  to  the  product  of  the  conductor 
length  and  resistance  per  unit  length.  The  resistance  per  unit  length  is  a  function  of  the 
conductor  material  and  conductor  cross-section. 

The  power  supply  power  is: 


Ppt  —  I^ind  +  (C.2) 

We  use  this  equation  in  chapter  five  to  characterize  the  power  needed  for  an  example 
application.  Sources  for  the  possibly  large  amount  of  vehicle  power  needed  include  solar 
photovoltaic  arrays,  nuclear  reactors,  and  earth-based  lasers. 
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Appendix  D.  Approximate  Mass  and  Moments  of  Inertia  for  the  POTV 


The  principal  moments  of  inertia  of  a  particular  POTV  design  (ASSET)  are  cal¬ 
culated  in  this  appendix.  The  resulting  expressions  are  used  in  the  MATLAB  controller 
design  and  performance  program  CONTROLLER.M. 

The  ASSET  POTV  system  with  Design  #2  conductor  configuration  consists  of  a 
modified  shuttle  external  tank  with  attached  rigid  conductor  towers.  We  model  the  system 
as  a  cylinder  with  attached  thin  rods.  The  rods  axe  assumed  to  have  constant  linear  mass 
density  p  (kg/m).  The  total  system  mass  is  then  the  simple  expression: 


Mt  —  Meyi  +  2p(/i  /2  +  A  -f-  Ici  +  ht) 
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Where  we  have  neglected  the  mass  of  the  guidewire  webbing,  PMGs,  power  generation 
equipment,  etc.  From  the  ASSET  study  (4)  we  find  that  the  mass  of  the  central  cylinder  is 
approximately  31,300  kg  (mass  of  an  empty  shuttle  external  fuel  tank).  By  no  means  are 
the  conductor  towers  required  to  all  have  the  same  linear  mass  density.  Indeed,  one  or  more 
conductor  towers  also  serving  as  docking  facilities  could  have  larger  linear  mass  densities 
than  the  other  towers.  Also,  towers  serving  as  structural  “anchors”  for  the  guidewire 
webbing  could  have  large  linear  mass  densities.  Here  we  assume  that  all  of  the  towers  have 
the  same  linear  mass  density  for  ease  of  derivation. 

The  approximate  moment  of  inertia  matrix  for  the  POTV  about  its  center  of  mass  is 
found  by  repeatedly  applying  the  parallel-axis  theorem  to  the  moment  of  inertia  matrices 
of  the  thin  rods.  The  parallel-axis  theorem  is  explained  by  Wiesel  (24:106-107).  The 
principal  moments  of  inertia  for  a  cylinder  and  a  thin  rod  are  found  in  the  dynamics  text 
by  Likins  (10:522-524).  Let  a  be  the  cylinder  radius  and  h  be  the  length  of  the  cylinder. 
After  simple  but  extended  calculation  the  moments  of  inertia  of  the  POTV  with  the  Design 
#2  conductor  configuration  are: 


A 

B 


-|-  Ip 


M,y,(3o"  +  h}) 


12 


Flat  (0.28  kg/m) 
I-beam  (0.76  kg/m) 


m 


Front  View: 


Figure  D.l.  Truss  Assembly  for  Conductor  Towers 
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End  and  side  views  of  the  proposed  truss  work  for  the  conductor  towers  in  shown 
in  figure  D.l.  We  assume  that  the  cross-sectional  area  of  each  tower  is  a  constsmt  one 
meter*.  Each  toww  section  consists  of  four  aluminum  “I-beams”  smd  twdve “flats”.  This 
material  is  readily  available  from  a  salvaged  external  tank  (4).  The  linear  mass  density  of  a 
single  I-beam  is  0.76  kg/m.  If  we  assume  that  the  “strip”  pieces  mentioned  in  the  ASSET 
report  can  be  divided  Imigthwise  into  eight  pieces,  then  the  linear  mass  density  of  the 
resulting  “flats”  is  0.28  kg/m.  With  these  assumptions  we  calculate  the  calculate  the  total 
linear  mass  density  oi  the  truss  to  be  approximately  7.3  kg/m.  This  is  the  figure  used  in 
the  CONTROLLER.M  program  listed  in  appendix  F.  A  smaller  total  linear  mass  density 
could  be  realized  if  the  tower  cross-sectional  area  specification  ( 1  meter*)  was  relaxed. 
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Here  is  the  dimension  data  for  the  design  used  in  the  chapter  five  examples: 


Mcyi  =  31300  kg 
p  =  7.3  kg/m 

11  =  310  m 

12  =  50  m 
It  =  50  m 

let  =  20  m 
le,  =  20  m 

The  resulting  total  system  mass  and  moments  of  inertia  are: 

Mt  =  37870  kg 
A  =  1.5374e  +  08  kg  •  m^ 
B  =  1.0236e  +  07kg-m^ 
C  =  1.5872c +  08  kg- m^ 

The  added  mass  due  to  the  conductor  towers  is  6670  kg. 
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Appendix  E.  Controller  Design  Program  Results 


E.  1  Steady-State  Tracking  Example 
»  cont 

WelcoBe  to  th«  POTV  Controllor  Dtsign  wd  Porfomeoica  Prograa. 

Altitude  of  reference  orbit  (kiloaeters)  ?  400 
Trig  angle  (u-tbeteg-laabda)  to  work  with  (degrees)  ?  33 
Mass  of  ASSET  excluding  tovers  (kg)  ?  31300 
Linear  aass  density  of  conductor  towers  (kg/a)  ?  7.3 
Length  of  xertical  (z)  in-plane  conductors  (aeters)  ?  50 
Length  of  horizontal  (y)  in-plane  conductors  (aeters)  ?  310 
Length  of  nonconducting  cross  support  towers  (aeters)  ?  50 
Length  of  vertical  (x)  cross  conductors  (aeters)  ?  20 
Length  of  horizontal  (y)  cross  conductors  (aeters)  ?  20 

Orbital  aean  notion  (rad/sec)  « 

1 . 13160-03 

Total  systea  aass  (kg)  » 

37870 

Principle  aonent  of  inertia  A  (kg-a~2)  “ 

1.5374e+08 

Principle  aoaent  of  inertia  B  (kg-B~2)  « 

1.0236e+07 

Principle  aoaent  of  inertia  C  (kg-B''2)  ^ 

1.5872e+08 

The  (FsySfGsys)  systM  is  cos^letely  controllable. 

The  (Fsys.Hsys)  systea  is  C(Mq>letely  observable. 
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Weighting  on  the  state  weighting  matrix  Q  ?  1 

Weighting  on  the  control  weighting  matrix  R  ?  1 

What  is  the  value  of  x  to  attain  and  track  (m)  ?  100 

What  is  the  value  of  y  to  attain  and  track  (m)  ?  100 

Wiat  is  the  value  of  z  to  attain  and  track  (m)  ?  100 

What  is  the  value  of  psil  to  attain  and  track  (rad)  ?  0 

What  is  the  value  of  psi2  to  attain  and  track  (rad)  ?  0 

What  is  the  value  of  psi3  to  attain  and  track  (rad)  ?  0 

What  is  the  initial  value  of  x  (m)  ?  100 

Wiat  is  the  initial  value  of  y  (m)  ?  100 

Vhat  is  the  initial  value  of  z  (m)  ?  100 

What  is  the  initial  veilue  of  psil  (rad)  ?  0 

What  is  the  initial  value  of  psi2  (rad)  ?  0 

What  is  the  initial  value  of  psi3  (rad)  ?  0 

How  long  of  a  time  response  to  plot  ?  10000 

Number  of  time  increments  (data  points)  ?  SOO 
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meters  meters 


Vertical  In-Plane  Position  (x)  vs  Time 


meters 


radians 


xlO-i*  Roll  Angle  (psi2)  vs  Time 


xlO-i* 
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Pitch  Angle  (psi3)  vs  Time 
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E.2  R-BAR  Docking  Approach  Example 
>>  cont 

Welcome  to  the  POTV  Controller  Design  and  Performance  Program. 

Altitude  of  reference  orbit  (kil(»aeters)  ?  400 
Trig  angle  (u-thetag-lambda)  to  work  with  (degrees)  ?  33 
Mass  of  ASSET  excluding  towers  (kg)  ?  31300 
Linear  mass  density  of  conductor  towers  (kg/m)  ?  7.3 
Length  of  vertical  (z)  in-plane  conductors  (meters)  ?  50 
Length  of  horizontal  (y)  in-plane  conductors  (meters)  ?  310 
Length  of  nonconducting  cross  support  towers  (meters)  ?  50 
Length*  of  vertical  (z)  cross  conductors  (meters)  ?  20 
Length  of  horizontal  (y)  cross  conductors  (meters)  ?  20 

Orbital  mean  notion  (rad/sec)  « 

1 . 1316e-03 

Total  systoB  mass  (kg)  > 

37870 

Principle  moment  of  inertia  A  (kg-n~2)  = 

1.5374e+08 

Principle  moment  of  inertia  B  (kg-m~2)  = 

1.0236e+07 

Principle  moment  of  inertia  C  (kg-m"2)  = 

1.5872e+08 


The  (Fsys,Gsys)  systm  is 
The  (Fsys.Hsys)  system  is 

Weighting  on  the  state  W' 
Weighting  on  the  control 


completely  controllable 
completely  observable. 

lighting  matrix  Q  ?  1000 
weighting  matrix  R  ?  1 


E-14 


What  is  the  value  of  x  to  attain  and  track  (m)  ?  50 

What  is  the  value  of  y  to  attain  and  track  (■)  ?  0 

What  is  the  value  of  z  to  attain  and  track  (■)  ?  0 

What  is  the  value  of  psil  to  attain  and  track  (rad)  ?  0 

What  is  the  value  of  psi2  to  attain  and  track  (rad)  ?  0 

What  is  the  value  of  psi3  to  attain  and  track  (rad)  ?  0 

What  is  the  initial  value  of  z  (■)  ?  100 

What  is  the  initial  value  of  y  (■)  ?  0 

What  is  the  initial  value  of  z  (■)  ?  0 

What  is  the  initial  value  of  psil  (rad)  ?  0 

What  is  the  initial  value  of  p8i2  (rad)  ?  0 

What  is  the  initial  value  of  psiS  (rad)  ?  0 

How  long  of  a  time  response  to  plot  ?  10000 

Number  of  time  increments  (data  points)  ?  500 
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meters 


xlO-3  Out-of-Plane  Position  (z)  vs  Time 


seconds 


xlO-*  Yaw  Angle  (psil)  vs  Time 
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E.3  V-BAR  Docking  Approach  Example 
»  cont 

Walcc»ie  to  the  POTV  Controller  Design  and  Performance  Program. 

Altitude  of  reference  orbit  (kilometers)  ?  400 
Trig  angle  (u-thetag-lambda)  to  work  with  (degrees)  ?  33 
Hass  of  ASSET  excluding  towers  (kg)  ?  31300 
Linear  mass  density  of  conductor  towers  (kg/m)  ?  7.3 
Length  of  wertical  (x)  in-plane  conductors  (meters)  ?  310 
Length  of  horizontal  (y)  in-plane  conductors  (meters)  ?  SO 
Length  of  nonccHiducting  cross  support  towers  (maters)  ?  50 
Length  of  vertical  (x)  cross  conductors  (meters)  ?  20 
Length  of  horizontal  (y)  cross  conductors  (meters)  ?  20 

Orbital  mean  motion  (rad/s«c)  = 

1.1316e-03 

Total  system  mass  (kg)  * 

37870 

Principle  moment  of  inertia  A  (kg-m''2)  « 

3.6521e+06 

Principle  moment  of  inertia  B  (kg-m"2)  • 

1.8873a+08 

Principle  moment  of  inertia  C  (kg-m''2)  « 

1.8702e+08 


Ihe  (FsySiGsys)  system  is 
The  (Fsys,Hsys)  system  is 

Veij^ting  on  the  state  wi 
Weighting  on  the  control 


co^}letely  cfmtrollable. 
coi^letely  observable. 

lighting  matrix  Q  ?  100 
weighting  matrix  R  ?  1 
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What  is  the  value  of  x  to  attain  and  track  (m)  ?  0 

What  is  the  value  of  y  to  attain  and  track  (m)  ?  50 

What  is  the  value  of  z  to  attain  and  track  (m)  ?  0 

What  is  the  value  of  psil  to  attain  and  track  (rad)  ?  0 

What  is  the  value  of  p3i2  to  attain  and  track  (rad)  ?  0 

What  is  the  value  of  psi3  to  attain  and  track  (rad)  ?  0 

What  is  the  initial  value  of  x  (m)  ?  0 

What  is  the  initial  value  of  y  (m)  ?  100 

What  is  the  initial  value  of  z  (m)  ?  0 

What  is  the  initial  value  of  psil  (rad)  ?  0 

What  is  the  initial  value  of  psi2  (rad)  ?  0 

What  is  the  initial  value  of  psi3  (rad)  ?  0 

Hov  long  of  a  tine  response  to  plot  ?  10000 

Miaber  of  tine  increnents  (data  points)  ?  500 
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Appendix  F.  Controller  Design  Program  Listing 


disp(’  ’); 

dispC ’Welcome  to  the  POTV  Controller  Design  end  Performance  Program.*); 
dispC’  ’); 
format  short  e 

Xdefine  constants 

G=6.67259e-ll;  ^gravitational  constant 

Me-5.976e24;  Xmass  of  the  earth  in  kg 

res6.378a06;  Xmean  equatorial  radius  of  earth  in  meters 

m«8.1el5;  Xmagnetic  dipole  moment  of  the  earth  in  Tesla 

d=-78.5*pi/180;  Xlatitude  of  austral  magnetic  pole  in  radians 

h>46.88;  Xlength  of  ET  modeled  as  a  cylinder  in  meters 
as4.206;  '/.radius  of  ET  modeled  as  a  cylinder  in  meters 

Xuser  enters  tower  lengths 

altitude=input(’ Altitude  of  reference  orbit  (kilometers)  ?  *); 
w*input(’Trig  angle  (u-thetag-lambda)  to  work  with  (degrees)  ?  ’); 
Metsinput(’Mass  of  ASSET  excluding  towers  (kg)  ?  ’); 
psinput(’Lineeu:  mass  density  of  c<»iductor  towers  (kg/m)  ?  ’); 
ll^inputC Length  of  vertical  (x)  in-plane  conductors  (meters)  ?  ’); 
12’:input(’ Length  of  horizontal  (y)  in-plane  conductors  (meters)  ?  ’); 
lt=input(’ Length  of  nonconducti^  cross  support  towers  (meters)  ?  '); 
lclsinput( ’Length  of  vertical  (x)  cross  conductors  (meters)  ?  ’); 
lc2»input(* Length  of  horizontal  (y)  cross  conductors  (meters)  ?  ’); 

Xunit  conversions 

r0«re4' altitude* 1000;  Xorbital  radius  of  reference  orbit  in  meters 
w*w*pi/180;  Xconvert  trig  angle  to  radians 
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Xcalculations 

n=sqrt(G*Me/rO“3) ;  Xorbital  nean  notion  (reference  orbit) 
dispC’Orbital  mean  motion  (rad/sac)  =*); 
disp(n) ; 

Mtot=Met+2*p*(ll+12+lcl+lc2+lt) ;  ‘/.Hass  of  entire  POTV  system  in  kg 
dispC’Total  systwa  mass  (kg)  =’); 
disp(Mtot) ; 

Iet=Met*([(a“2/2)  0  0;0  ((3*a‘2+h“2)/12)  0;0  0  ((3*a“2+h“2)/12)] )  ; 
Ill=2*ll*p*([0  0  0;0  ((ll‘2/12)+((ll+h)/2)“2)  0;0  0 

((ir2/12)  +  ((ll+h)/2)“2)]); 

I12=2*12*p*([((12-2/12)+((12/2)+a)“2)  0  0;0  0  0;0  0 
((12‘2/12)+((12/2)+a)“2)]); 

Ilt=2*lt*p*(C((lt“2/12)+((lt/2)+a)“2)  0  0;0  ((lt“2/12)+((lt/2)+a)“2)  0; 

0  0  0]); 

Ilcl=«2*lcl*p*([((lt+a)“2)  0  0;0  ((lcl'-2/12)+(lt+a)“2)  0;0  0  (lcl“2/12)] ) ; 
Ilc2=2*lc2*p*([((lc2''2/12)+(lt+a)‘2)  0  0;0  (lt+a)'2  0;0  0  (lc2“2/12)]) ; 
I=Iet+Ill+I12+Ilt+Ilcl+Ilc2;  XHonent  of  inertia  matrix  in  kg*B"2 
A=I(1,1);  Xprinciple  moments  of  inertia 
B=I(2,2); 

C*I(3,3); 

disp(’ Principle  moment  of  inertia  A  (kg-m'‘2)  =’); 
di8p(A); 

di8p( 'Principle  moment  of  inertia  B  (kg-n"2)  *’); 
di8p(B) ; 

di8p( 'Principle  moaMnt  of  inertia  C  (kg-m''2)  -’); 
di8p(C); 

Xeet  up  8y8tem  matricee  F,G,H,J  (open-loop) 

XsyatMi  matrix 
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Fsys=[0  10000000000; 

(3*n''2)  0  0  (2*n)  0000000  0; 

000100000000; 

0  (-2*n)  0000000000; 

000001000000; 

0000  (-(n“2))  0000000; 

00000001000  0; 

000000  (-(n“2)*(C-B)/A)  0  0  (-n*(-A-B+C)/A)  0  0; 
000000000100; 

0000  ((3*n*2)*(A-C)/(rO*B))  0  0  (-n*(A+B-C)/B)  ((4*n“2)*(A-C)/B)  000; 
00000000000  1; 

0  (3*n/(2*r0))  ((3*n'-2)*(B-A)/(rO*C))  0  0  0  0  0  0  0  ((-3*n“2)*(B-A)/C)  0]  ; 
Xinput  aatrix 

G8ys®(l/r0“3)*([0  00000000; 

0  0  (-■♦12*siB(d)/Mtot)  (-m*12*8in(d)/Mtot)  0  0  (-■♦lc2*8in(d)/Htot) 
(-■*lc2*8in(d)/Mtot)  0; 

00000000  0; 

(■*ll*»in(d)/Mtot)  (B*ll*sin(d)/Htot)  0  0  (B*lcl*sin(d)/Mtot) 
(■*lcla8iB(d)/Ntot)  000; 

00000000  0; 

(■<‘11*008  (d) *8111(8) /Mtot)  (■*ll*co8(d)*8in(v)/Ntot) 

(-2*b*12*co8 (d) *cos (w) /Mtot )  ( -2*m*12*co8 (d) *co8 (w) /Mtot ) 
(■*lci*co8(d)*siii(w)/Mtot)  (n*lcl*co8(d)*8iB(v)/Mtot) 

(-2*B*lc2*cos(d)*co8(B)/Mtot)  (-2*B*lc2*co8(d)*co8(w)/Mtot)  (rO“3/Mtot) ; 
00000000  0; 

0  0  (-■*12“2*co8(d)*co8(H)/A)  (■*12"2*co8(d}*co8(w)/A) 
(-■*lt*lcl*8in(d)/A)  (B*lt*lcl*8in(d)/A)  000; 

00000000  0; 

(-■*ll“2*co8(d)*8iB(w)/(2*B))  (■*ll''2*co8(d)*8in(w)/(2*B))  0000 
(-■*lt*lc2*8in(d)/B)  (■*lt*lc2*8in(d)/B)  0; 
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000000000; 

(m*ll*2*sin(d)/C2*C))  (-■*ll“2*sin(d)/(2*C))  (■*12“2*sin(d)/(2*C)) 
(-ni*12“2*sin(d)/(2*C))  0  0  0  0  0]); 

^output  matrix  (selects  positions  and  attitude  angles  only) 

Hsys=[l  00000000000; 

001000000000; 

000010000000; 

000000100000; 

000000001000; 

00000000001  0]; 


Jsys=zeros(6,9) ;  Xno  fed-through  control  signals 

Xtest  controllability  and  observability 
disp(’  ’); 

Hc-ctrb(Fsys,Gsy8) ; 
if  rank(Mc)*=12, 

dispC’The  (Fsys.Gsys)  systMi  is  completely  controllable.’): 

else 

disp(’The  (Fsys.Gsys)  system  has  uncontrollable  modes!’); 

end 

disp(’  ’); 

Mo>ob8v(Fsys .Hays) ; 
if  rank(Mo)«»12  , 

disp(’'I%e  (Fsys.Hsys)  system  is  c^pletely  observable.’); 

else 

disp(’The  (Fsys.Hsys)  systMi  has  unobservable  modes!’); 

end 


F-4 


Xderiva  LQR  coapansator  K 

Xuser  enters  state  and  control  signal  weightings 


Xnormalized  state  weighting  matrix 

q« input (’Weighting  on  the  state  weighting  matrix  Q  ?  ’); 
Q=q*([l  00000000000; 

010000000000; 

001000000000; 

000100000000; 

000010000000; 

000001000000; 

000000  rO  00000; 

OOOOOOOrOOOOO; 

00000000  rO  000; 

000000000  rO  00; 

0000000000  rOO; 

00000000000  rO]); 

r>input( ’Weighting  on  the  control  weighting  matrix  R  ?  ’); 
Rnr’t'eyeO) ;  Xequal  weighting  on  each  control 

K>lqr(Fs3rB,6s]rs,q,R) ; 

Xderiwe  non-square  pre-filter  N  using  optimal  pseudo-inverse 
N«-pinv (Hsys* inv (Fsjs-GsyseK) *Gsys ) ; 

Xsetup  closed- loop  system  matrices 
Fcl^Fsys-GsyseK ; 

Gcl>Gsys*M; 

Hcl«llsys ; 

Jcl^EerosCG) ; 
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^linear  simulation  of  states  due  to  initial  condition  and  reference  cmd 
rx^inputC’What  is  the  value  of  z  to  attain  and  track  (m)  ?  ’): 

ry»input(’What  is  the  value  of  y  to  attain  and  track  (m)  ?  ’); 

rz-input( 'What  is  the  value  of  z  to  attain  and  track  (m)  ?  ’); 

rpsil=input( 'What  is  the  value  of  psil  to  attain  zuid  track  (rad)  ?  ’); 

rpsi2=input( 'What  is  the  value  of  p8i2  to  attain  zuid  track  (rad)  ?  ’); 

rp8i3-input( 'What  is  the  value  of  psiS  to  attain  and  track  (rad)  ?  ’); 

xO* input (’What  is  the  initial  value  of  x  (■)  ?  ’); 

y0= input ( 'What  is  the  initial  value  of  y  (m)  ?  ’); 

zO=input( 'What  is  the  initial  value  of  z  (m)  ?  ’); 

psilO=input(’What  is  the  initial  value  of  psil  (rad)  ?  ’); 

psi20=input('What  is  the  initial  value  of  psi2  (rad)  ?  ’); 

p8i30«:input(’What  is  the  initial  value  of  p8i3  (rad)  ?  ’); 

X0*[x0  0  yO  0  zO  0  psilO  0  psi20  0  psiSO  0]’; 
tmax«input ( ’ How  long  of  a  time  response  to  plot  ?  '); 
ntinc^B  input  (’Number  of  time  increments  (data  points)  ?  ’); 
t®0 : (tmax/ntinc) : tmax ; 
tleng=length(t) ; 

r=zeros(tleng,6) ;  Xconstruct  reference  command 

r(: ,l)srx*ones(tleng,l) ; 

r(: ,2)sry*ones(tleng,l) ; 

r(: ,3)»rzeones(tleng,l) ; 

r(: ,4)srpsiieone8(tleng,l); 

r(:  ,5)>rpsi2'*ones(tleng,l) ; 

r( : ,6)Brp8i3eoBe8(tleng,l) ; 

Cy,x]>lsim(Fcl,Gcl,Hcl, Jcl,r,t,X0) ;  Xoutpats  and  states 
u*(M*r’-K*x’) ' ;  Xcontrol  vector  includes  8  currents  and  thrusters 
f «(Gsyseu’ ) ’ ;  Xforces 
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plot(t,y(: ,1)) ; 


title( ’Vertical  In-Plane  Position  (x)  vs  Time’); 

xlabel ( ’ seconds ’ ) ; 

ylabel ( ’meters ’ ) ; 

meta  x; 

pause; 

plot (t , (Mtot*f ( : , 2) ) ) ; 

titleC’Force  in  Vertical  Direction  (x)  vs  Time’); 

xlabel ( ’ seconds ’ ) ; 

ylabel ( ’ Newtons ’ ) ; 

meta  fx; 

pause ; 

plot(t,y(: ,2)) ; 

title (’Horizontal  In-Plane  Position  (y)  vs  Time’); 

xlabel ( ’ seconds ’ ) ; 

ylabel (’meters’) ; 

meta  y; 

pause ; 

plot(t , (Mtot*f ( : ,4))) ; 

title(’Force  in  Horizontal  (y)  Direction  vs  Time’); 

xlabel ( ’ seconds ’ ) ; 

ylabel (’Newtons’) ; 

meta  fy; 

pause ; 

plot(t,y(: ,3)) ; 

title(’Out-of -Plane  Position  (z)  vs  Time’); 
xlabel ( ’ seconds ’ ) ; 
ylabel ( ’meters ’ ) ; 
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meta  z; 


pausa : 

plot(t , (Mtot*f ( : ,6) ) ) ; 

titla( ’Force  in  Out-of -Plane  (z)  Direction  ra  Tine’) 
zlabel ( ’ seconds ’ ) ; 
ylabel ( ’ Heatons ’ ) ; 

■eta  fz; 

pause ; 

plot(t,y(: ,4)) ; 

titleC’Yaw  Angle  (psil)  vs  Time’); 
zlabel ( ’ seconds ’ ) ; 
ylabel ( ’ radians ’ ) ; 

■eta  psil; 

pause ; 

plot(t,(A*f (: ,8))); 
titleC’Moaent  About  z-Azis  vs  Tine’); 
zlabel ( ’ seconds ’ ) ; 
ylabel (’■-■’); 

■eta  ■z; 
pause ; 

plot(t,y(: ,5)) ; 

title(*Roll  Angle  (psi2)  vs  Tine’); 
zlabel ( ’ seconds ’ ) ; 
ylabel ( ’ radians ’ } ; 

■eta  pBi2; 
pause : 
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plot(t, (B*f ( : ,10))) ; 

titleC’Moraent  About  y-Axis  ts  Tine’); 

xlabel ( ’ seconds ’ ) ; 

ylabel(’N-m’); 

meta  my; 

pause : 

plot(t,y(: ,6)); 

title( ’Pitch  Angle  (psiS)  vs  Tine’); 

xlabel ( ’ seconds ’ ) ; 

ylabelC ’radians’) ; 

meta  psi3; 

pause ; 

plot(t,(C*f (: ,12))) ; 

titleC ’Moment  About  z-Axis  vs  Tine’); 

xlabeK ’seconds’ ) ; 

ylabeK’M-n’) ; 

meta  mz; 

pause ; 

plot(t,u(: ,1)); 

title (’Current  Control  Signal  (Ilu)  vs  Time’); 

xlabel ( ’ seconds ’ ) ; 

ylabel ( ’ amperes ’ ) ; 

meta  ilu; 

pause ; 

plot(t,u(: ,2)); 

titleC ’Current  Control  Signal  (lid)  vs  Time’); 
xlabel ( ’ seconds ’ ) ; 
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plot(t,u(: ,7)) : 

title (’Current  Control  Signal  (I2n)  vs  Time’); 
zlabel ( ’ seconds ’ ) ; 
ylabel ( ’ amperes ’ ) ; 

■eta  i2n; 
pause ; 

plot(t,u(: ,8)) ; 

title (’Current  Control  Signal  (128)  vs  Time’); 

xlabel ( ’ seconds ’ ) ; 

ylabel ( ’ amperes ’ ) ; 

meta  i2s; 

pause ; 

plot(t,u(: ,9)) ; 

title (’Thruster  Control  Signal  (T)  vs  Tijie’); 
xlabel ( ’ seconds ’ ) ; 
ylabel ( ’ newtons ’ ) ; 
meta  t; 
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ylab«l ( ’ aBp«r«3 ’ ) ; 
meta  ild; 
pause ; 

plot(t,u(: ,3)) ; 

title (’Current  Control  Signal  (I2e)  vs  Time’); 

xlabeK  ’  seconds  ’ )  ; 

ylabel ( ’ amperes ’ ) ; 

meta  i2a: 

pause ; 

plot(t,u(: ,4)) ; 

title  (’Current  Control  Signal  (I2ir)  vs  Time’); 

xlabel ( ’ seconds ’ ) ; 

ylabel ( ’ amperes ’ ) ; 

meta  i2«; 

pause ; 

plot(t,u(: ,5)) ; 

title (’Current  Control  Signal  (I In)  vs  Time’); 
xlabel ( ’ seconds ’ ) ; 
ylabel  ( ’  wperes  ’ )  ; 
meta  iin; 

pause ; 

plot(t,u(: ,6)); 

title (’Current  Control  Signal  (Ils)  vs  Time’); 

xlabel ( ’ seconds ’ ) ; 

ylabel ( ’ amperes ’ ) ; 

meta  ils; 

pause ; 
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